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1 Introduction 

In this Chapter, we quickly review some basic definitions and concepts from thermodynamics. We 
then provide a brief description of the first and second laws of thermodynamics. Next, we discuss the 
mathematical consequences of these laws and cover some relevant theorems in multivariate calculus. 
Finally, free energies and their importance are introduced. 


1.1 Basic concepts 

1.1.1 State function versus path function 

A state function is a function that depends only on the current properties of the system and not on the 
history of the system. Examples of state functions include density, temperature, and pressure. 

A path function is a function that depends on the history of the system. Examples of path functions 
include work and heat. 


1.1.2 Intensive property versus extensive property 

An extensive property is a characteristic of a system that is proportional to the size of the system. 
That is, if we double the size of the system, then the value of an extensive property would also 
double. Examples of extensive properties include total volume, total mass, total internal energy, etc. 
Extensive properties will be underlined. For example, the total entropy of the system, which is an 
extensive property, will be denoted as S. 

An intensive property is a characteristic of a system that does not depend on the size of the system. 
That is, doubling the size of the system leave the value of an intensive property unchanged. Examples 
of intensive properties are pressure, temperature, density, molar volume, etc. By definition, an inten- 
sive property can only be a function of other intensive properties. It cannot be a function of properties 
that are extensive because it would then depend on the size of the system. 


1.2 Brief review of thermodynamics 
1.2.1 The first law of thermodynamics 

The first law of thermodynamics is simply a statement of the conservation of energy. Energy can take 
on a variety of forms, for example kinetic energy, chemical energy, or thermal energy. These different 
forms of energy can transform from one to another; however, the sum total of all the types of energy 
must remain constant. 
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Let’s apply the first law of thermodynamics to a closed system (i.e. a system that can exchange heat 
and work with its surroundings, but not matter). The first law for a closed system can be written as 


dU = SQ - SW H 


(u) 


where U is the internal energy of the system, SQ is the heat (thermal energy) transferred to the system, 
and SW is the work performed by the system. Other forms of energy may contribute to the energy 
balance, such as kinetic energy or potential energy (e.g., from a gravitational or electrostatic field). 

1.2.2 The second law of thermodynamics 

The second law of thermodynamics formalizes the observation that heat is spontaneously transferred 
only from higher temperatures to lower temperatures. From this observation, one can deduce the 
existence of a state function of a system: the entropy S. The second law of thermodynamics states 
that the entropy change dS of a closed, constant- volume system obeys the following inequality 



( 1 . 2 ) 


where T is the absolute temperature of the system, and SQ is the amount of heat transfered to the 


system. A process will occur spontaneously in a closed, constant- volume system only if Eq. (1.2) 
is satisfied. For a reversible process, the equality is satisfied; for an irreversible process, the entropy 
change is greater than the right-hand side of the relation. 

Note that the second law of thermodynamics is unique among the various laws of nature in that it 
is not symmetric in time. It sets a direction in time, and consequently there is a distinction between 
running forward in time and running backwards in time. We can notice that a film is being played in 
reverse because we observe events that seem to violate the second law. 

1.3 The fundamental equation of thermodynamics 

Now consider a closed system that can alter its volume V '. In this case, the work performed by the 
system is SW = pdV . Combining the first and the second laws of thermodynamics for a closed 
system (i.e. inserting the inequality in Eq. (1.2) into Eq. (1.1)), we obtain 


dU < TdS — pdY for constant N 


(1.3) 


For any spontaneous change (process) in the system, the inequality given in Eq. (1.3) will be satisfied. 
The equality will be satisfied only in a reversible process. 

An isolated system is a system that does not exchange work SW = 0, heat SQ = 0, or matter 
dN = 0 with its surroundings. Consequently, the total internal energy and volume remain constant; 
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that implies that dU = 0 and dV = 0. Substituting these relations into Eq. (1.3), we find that 
processes occur spontaneously in an isolated system only if the entropy does not decrease. In this 
case, 

dS> 0 (1.4) 

Note that in an isolated system, every spontaneous event that occurs always increases the total entropy. 
Therefore, at equilibrium, where the properties of a system no longer change, the entropy of the 
system will be maximized. 

For a system where entropy and volume are held fixed (i.e. dS = 0 and V = 0), a process will occur 
spontaneously if 

dU < 0 at constant S, Y, and N (1.5) 

For a reversible process, where the system is always infinitesmally close to equilibrium, the equality in 
Eq. (1.3) is satisfied. The resulting equation is known as the fundamental equation of thermodynamics 

dU = TdS — pdY at constant N (1.6) 
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1.4 The calculus of thermodynamics 

From the fundamental equation of thermodynamics, we can deduce relations between the various 
properties of the system. To see this, let’s consider a function / with independent variables x and y. 
The differential of / (i.e. the total change in /) can be written as: 




dx+ (%) d y 


(1.7) 


The first term represents the change in / due to changes in the independent variable x, and the second 
term representce changes due to the independent variable y. Note that Eq. (1.7) is just a generalization 
of a first order Taylor series expansion to a function of two variables. 


If we consider the internal energy of the system U to be a function of the variables S and V, then 
taking / — ► U , x — ► S, and y — > V, we find 


dU = 


dU 

dSsv 


dS+l §) s dV 


( 1 . 8 ) 


Comparing Eq. (1.8) with the fundamental equation of thermodynamics Eq. (1.6), we can make the 
identifications 


*L) = T 

dSjy 

dU' 
dV 


~P 


(1.9) 

( 1 . 10 ) 


Therefore, we see that the temperature and pressure of the system can be related to derivatives of its 
internal energy. 

For most functions, the order of differentiation does not matter. That is 


d_(d£ 

dx \dy 


xJ y 


d_fd£\ 

d V \ dx ) y 


If we apply Eq. (1.1 1) to the internal energy, we find 


9 1 

( dU N 

\ 


' d 

(dU\ ' 

dv 1 


K_ 

s 

dS 1 

\9YJs_ 


dT\ 
WJs 


dp 

dS 


v 


( 1 . 11 ) 


( 1 . 12 ) 


where we have used Eqs. (1.9) and (1.10). These types of relations are known as Maxwell relations. 
We will encounter more of these kind of relations later on. 


There are three additional relations that need to be mentioned. These relations are useful in converting 
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properties that depend on “unmeasureable” quantities, such as entropy, to properties that are measure- 
able, such as temperature or pressure. The first is a generalization of the chain rule to functions of 
multiple variables 


d£ 

dy 


0 £ 

dz 


dz 

dy 


(1.13) 


To determine the other two relations, let’s consider a function of three variables that is constrained to 
be equal to zero. That is 

f(x,y,z ) = 0. 

This defines a two-dimensional surface embedded in a three-dimensional space. The above equation 
can be interpreted as defining the functions: 

x = x(y, z) y = y(x, z) z = z(x, y). 


Each of these functions can be expanded in terms of its respective independent variables 


dx 

dy 

dz 



dy + 

Z 

dx + 

z 

dx + 

v 



(1.14) 

(1.15) 

(1.16) 


By substituting Eq. (1.15) in to Eq. (1.14) to eliminate the dy term, we find 


dx = 


0 = 



dy_ 

dx 


dx + 


dy 

dz 


dz 


+ 




dx + 




(1.17) 

(1.18) 


This equation should hold for any value of dx and dz. In order for this to be true, the coefficients of 
the dx and dz must vanish. For the dx coefficient, we find 


fdx\ = (dy \~ 1 
\0yJ z \dx) z 

The coefficient of the dz term leads to 

f dx\ f dx\ f dy\ 

\d~z) y = ~\'d~y) z \d~z) 

This relation is known as the triple product rule. 


(1.19) 


( 1 . 20 ) 
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1.5 Open systems 


We can extend Eq. (1.6) to open systems (i.e. systems in which the number of moles N can vary) by 
including a term called the chemical potential fi. Mathematically, this quantity represents the increase 
in the internal energy when a small amount of material is introduced to the system at constant total 
entropy and volume: 


fi = 



( 1 - 21 ) 


The fundamental equation of thermodynamics can then be written as 


dU = TdS - pdY + fidN 


( 1 . 22 ) 


Let’s determine what the chemical potential is. To do this, we rewrite all total quantities in terms of 
molar properties. For example, the total internal energy can be written as U = NU , where U is the 
molar internal energy of the system, and N is the total number of moles in the system. Substituting 
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these relations into the fundamental equation, we find 

d(NU) = Td(NS) - pd(NV) + pdN 
NdU + UdN = NTdS + TSdN - NpdV - pVdN + pdN 
NdU = N(TdS - pdV) - (U +pV -TS - p)dN 

dU = TdS - pdV - L (G - p)dN (1 .23) 

where in the last line, we have introduced the definition of the molar Gibbs free energy G = U + 
pV - TS. 


The molar internal energy U is an intensive property of the system; therefore, it should be independent 
of extensive properties of the system, in particular, the total number of moles in the system N. In 
order for this to be true, the chemical potential /i must be equal to the molar Gibbs free energy G. In 
otherwords: 

p = G (1.24) 

Note that this derivation is restricted to pure substances. For multicomponent systems, we need to 
generalize this relation. This will be done later. 


1.6 Legendre transforms and free energies 

The natural variables of the internal energy U are the entropy S, volume V, and total number of moles 
N of the system. In many situations, however, these variables are not convenient. 

We can easily arrive at a new function that has different natural variables by performing a Legendre 
transform. For example, to arrive at a new state property that posesses the independent variables T 
and V, we define the Helmholtz free energy A as: 

A = U -TS (1.25) 

Inserting this relation into Eq. (1.3), generalized to open systems, we find 

dU < TdS - pdV + pdN 
d(A + TS) < TdS - pdV + pdN 

dA < —SdT - pdV + pdN (1.26) 

From this equation, we see that for a system with the temperature, volume, and total number of moles 
held fixed (i.e., dT = 0, dV = 0, and dN = 0), a process is spontaneous if it decreases the Helmholtz 
free energy. 
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In addition, at equilibrium where the equality holds, we find 

dA = —SdT - pdY + pdN (1.27) 

From this expression, we see that the natural variables of the Helmholtz free energy A are the tem- 
perature, volume, and total number of moles of the system. 

Similarly, if we define the Gibbs free energy G = U — TS + pY, then the fundamental equation of 
thermodynamics becomes 


dG = -SdT + Ydp + pdN (1.28) 

The Gibbs free energy is minimized for a system at constant temperature, pressure, and total number 
of moles. The Gibbs free energy is important because in most experiments the temperature and 
pressure are variables that we control. This will become useful to us later when we consider phase 
equilibria. 

The corresponding equation for the enthalpy H = U + pY is 

dH = TdS + Ydp + pdN (1.29) 


As we have seen, free energies such as the internal energy and Gibbs free energy are useful in that 
they tell us whether a process will occur spontaneously or not. A process in which the requisite free 
energy decreases will occur spontaneously. A process in which the free energy increases will not 
occur spontaneously. This does not mean that the process cannot happen; we can force the process 
to occur by performing work on the system. Therefore, we see that free energies are useful to us, 
qualitatively, in that they tell us the direction in which things will naturally happen. 

Free energies also provide us with quantitative information about processes. The change in the free 
energy is equal to the maximum work that can be extracted from a spontaneous process, or in the case 
of a non-spontaneous processes, the minimum amount of work that is required to cause the process 
to occur. 
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Free energies also have an additional, fundamental importance. Once the mathematical form of the 
free energy of a system is known in terms of its natural independent variables (e.g., the Gibbs free 
energy as a function of T, p, and N), then all the thermodynamic properties of the system can be 
determined. In the remainder of the course, we will be learning how to both develop approximate 
models for the free energy and how to use these models to estimate the thermodynamic behavior of 
various systems. 
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2 Single component systems 

In this Chapter, we describe the basic thermodynamic properties of single component systems. We 
begin with a qualitative description of their general phase behavior. Then, we discuss the mathemati- 
cal relations that govern this behavior. 


2.1 General phase behavior 

Figure 2. 1 shows the pressure-temperature projection of the phase diagram for a general one-component 
system. Depending on the temperature and pressure, the system can exist in either a solid, liquid, or 
vapor phase. Lines separate the various phases. On the lines, two phases coexist. The line separating 
the vapor and liquid phases is known as the vapor pressure curve. On crossing this curve, the system 
will transform discontinuously from a liquid to a vapor (or vice-versa). At high temperatures, the va- 
por pressure curve ends at a critical point. Beyond this point, there is no real distinction between the 
vapor and liquid phases. By going around the critical point, a liquid can be continuously transformed 
into a vapor. 

The line separating the solid and liquid phases is known as the melting or freezing curve. The line 
separating the solid and vapor phases is known as the sublimation curve. The point where the vapor 
pressure curve, the melting curve, and the sublimation curves meet is the triple point. At these 
conditions, the solid, liquid, and vapor phases can simultaneously coexist. 

In Figure 2.2, we show the temperature-density phase diagram for a general pure substance. As 
with the pressure-temperature diagram, the temperature-density phase diagram is divided by various 
curves into vapor, liquid, and solid phases. Outside these curves, the system exists as a single phase. 

Inside these curves, multiple phases coexist. For example, if a system is prepared at a state corre- 
sponding to point A in the diagram, it will divied into a vapor phase, with density p^°\ and liquid 
phase, with density p® . From the phase diagram, we can also determine the relative amounts of the 
coexisting phases. Let’s consider a system consisting of N total moles that are separated into a liquid 
phase that occupies a volume and a vapor phase that occupies a volume V^ 9 \ From a mole 
balance, we have 


N = P {1) ¥ {1) + p {a) v {9) 
p(¥ (l) + y(g)) = p V)y( 0 + p (g)y(g) 

(p {l) -p)¥ (l) = (p-p i9) )¥ i9) 

Y «) _ p - p te) 

¥® p®-p K ’ 

where p — N/ ( \ : ; is the overall density of the system. Equation (2.1) is known as the lever 

rule. It allows us to compute the relative amounts of the two-coexisting phases. 
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The dashed-line represents the triple point. Anywhere along the dashed-line, the vapor, liquid, and 
solid phases can simultaneously exist. 


2.2 Conditions for phase equilibrium 

Now let’s derive the mathematical conditions for equilibrium between two coexisting phases. We 
consider an isolated system that is separated into two phases, which we label A and B. The volume 
occupied by each phase can change; in addition, the both phases can freely exchange energy and 
material with each other. Because the system is isolated, the total energy U, the total volume Y, 
and the total number of moles N in the system must remain constant. This leads to the following 
relations: 


Z/( A) + J7 (B) = u 
y( A ) + y( B ) = y 
+ 7V (b) = N 


dU {A) = -dU {B) 
dV {A) = -dV {B) 
dN {A) = -dN^ B) 


( 2 . 2 ) 

(2.3) 

(2.4) 


where C/M is the total energy of phase i, Y^> is the total volume of phase i, and is the total 
number of moles in phase i. 
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The total entropy of an isolated system at equilibrium is maximized. Therefore, we have 

dS = 0 (2.5) 

The change in the total entropy of the system is given by the sum of the entropy change in phase A 
and the entropy change in phase B. This leads to 

dS = dSM + d§W 

i A A ) ,,A) 

= AjrdUW + ^-rjrdV^ A) ~ ^jjrdN^ 

t ( a ) _ r( A ) _ t( a ) 

1 n( B ) ,,(B ) 

+ ~^rdU {B) + - ^dN^ (2.6) 

T( b ) - t( b ) ■ T( B ) 

where we have rearranged Eq. (1.22) to obtain how the entropy of a phase changes with its energy, 
volume, and number of moles. 
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Substituting the relations given in Eq. (2.4) into Eq. (2.6), we find 


1 A A ) //( A ) 

dS = —r^dU {A) + ^jjTdV iA) - ^j-rrdNW 

T{A) - t(A) - t(a) 

ij^dUW - ^-rdV {A) + ^rrdN^ A) 

T(B) - rp(B) - J-(S) 


0 = 


1 


1 


J'(A) T(B) 


dU {A) + 


P 


(A) 


P 


(B) 


T(A) t ( b ) 


(A) 


d y(A) _ Mil _ 


M 


(B) 


T(A) t ( b ) 


dN (A ) (2.7) 


The quantities dU^ A \ dV ( A \ and dN (A) on the right-hand side of Eq. (2.7) can be chosen arbitrarily. 
From Eq. (2.6), we know that the left-hand side of Eq. (2.7) must equal zero. The only manner in 
which to guarantee the equality between the two sides of Eq. (2.7) is if coefficients of the dU^ A \ 
dV^ A \ and dN^ terms are each equal to zero. As a result, this implies the relations: 

rp(A) _ j:(B) 

P (A) = P (B) (2.8) 


Therefore, we find that the temperatures, pressures, and chemical potentials are equal for coexisting 
phases in equilibrium. We will later extend this derivation to multicomponent system. 


One interesting point to note is that the labels A and B which we used to derive Eqs. (2.8) do not 
have to refer to different phases. For example, A and B can refer to different parts of a one phase 
system. Therefore, Eqs. (2.8) can be interpreted as stating the the temperature, pressure, and chemical 
potential of a system at equilibrium are uniform (N.B., we did not include the influence of external 
fields, such as gravitional or electrostatic fields). 


2.3 The Clapeyron equation 

In this section, we derive the Clapeyron equation. This equation relates changes in the pressure to 
changes in the temperature along a two-phase coexistence curve (e.g., the vapor pressure curve or the 
melting curve). Note that the condition for equilibrium between two phases is given by 

^ ^ 

q{A) = q{B) 

dG {A) = dG^ B) 

-S iA) dT + V {A) dp = -S {B) dT + V {B) dp 

dp _ gw - sm 

dr vw - vm ’ 

This is one form of the Clapeyron equation. It relates the slope of the coexistence curve to the entropy 
change and volume change of the phase transition. 
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Entropy is not directly measureable, and, therefore, the Clapeyron equation as written above is not 
in a convenient form. However, we can relate entropy changes to enthalpy changes, which can be 
directly measured. At equilibrium, we have 

^) = ^) 
q{A) = g (B) 

_ ts ( A ) = H ^ — TS 

S( A ) - = i( H (A) - tf (B) ) (2.10) 

Thus, the entropy change of a phase transition, which is not directly measureable, can be determined 
from the enthalpy change of the phase transition, which is directly measureable. 


Substituting this relation into Eq. (2.9), we find 

dp _ - H 

dT ~ T(1/( A ) - l/W) 

This is the more commonly used form of the Clapeyron equation. 


( 2 . 11 ) 
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3 Multicomponent systems 

In this section, we examine the thermodynamics of systems which contain a mixture of species. 
First, we generalize the thermodynamic analysis of the previous section to multicomponent systems, 
deriving the Gibbs phase rule. Then we describe the general phase behavior of binary and ternary 
mixtures. 


3.1 Thermodynamics of multicomponent systems 

3.1.1 The fundamental equation of thermodynamics 


In this section, we extend the results of the previous lectures to multicomponent systems. All that 
needs to be done is to define a chemical potential for each species a in the system 


dU = TdS — pdY + p a dN a 

a 


(3.1) 


where p a is the chemical potential of component a. From this, we see that 


Ma = 


f— ) 

\ 9N a ) s,Y,N a i-£ a 


(3.2) 


Physically, p a is the change in the internal energy of the system with respect to an increase in the 
number of moles of species a, while holding all number of moles of all other species constant. 
The other forms of the fundamental equation of thermodynamics can be similarly generalized by 
performing the Legendre transform: 


dA = —SdT - pdY + VadN a 

a 

dG — — SdT -\- V dp T ^ ^ p^dN^ 

a 

dH = TdS + Ydp + p a dN a (3.3) 

a 


From these 
potential: 


relations, we see that there are the following alternate interpretations of the chemical 


Ma = 


\9N a J T ,y,N aljil 



T ->P ->N a t -£oc 





(3.4) 


3.1.2 Phase equilibria 

The conditions for phase equilibria can also be extended to multicomponent systems. Consider an 
isolated system consisting of uj components and two phases, which we refer to as A and B. The 
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system is isolated, with a total internal energy of U, a total volume of Y, and N a moles of species a. 
Because the system is isolated, we have the following 

y(A) + U (B) = jj ^ dU {A) = ~dU {B) 

y{A) + y(B) = y ^ d y{A) = _ d y( B ) (3.5) 

Ni A) +Ni B) =N a =► dN {A) =-dNi B) 

where U ij ) is the internal energy of phase i, is the volume of phase i, and A)/ is the number of 

moles of species a in phase i. 

For an isolated system at equilibrium, the total entropy is maximized. Stated mathematically, we have 

dS = dSM + d§W = 0 

1 n(A) ,X A ) 

= -h,duW + 4T dV (A) _ y 1^-dN^ 

y(A) ~ ~ ZL-/ r(A) a 

a 

OL 

Inserting the constraint relations given in Eq. (3.5) into Eq. (3.6), we find 

ds= ( A du (A) + ( P^_ _ dv (A) _ y ( _ 

\T ( a ) T( B ) ) ~ \T( a ) T( b ) J - 

This is equal to zero only if all the coefficients of the change terms are zero, 
find 


(3.6) 


(B) \ 

^JdNi A) (3.7) 
As a consequence, we 


y( A ) — y( B ) 

p (A)= p (B) 

fi ^ for all components a 

This argument can be generalized to a system containing n phases and uo components. In this case, 
we have the temperature, pressure, and chemical potentials of each species are equal in each phase. 

y(A) _ ji(B) — ... — j»(?r) 

p ( A )= p ( B ) = ...= p W 

fi^ = fj,^ = • • • = fi ^ for all components a (3.8) 


3.1.3 Gibbs phase rule 

How many variables need to be specified in order to fix the state of a system? In order to fix the 
state of a one-phase system, the composition of the phase must be specified as well as two additional 
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intensive variables (e.g., temperature and pressure). For a phase with lu components, cu — 1 mole 
fractions are required to specify the composition. Therefore, a total of uj + 1 intensive variables are 
required to specify the state of a single phase. 

For a system with n phases, there are a total of (u + l)7r unknowns. However, not all of these are 
independent. The conditions for phase equilibrium (see Eqs. (3.8)) give us (cu + 2)(tt — 1) equations 
that must be satisfied between each of the phases. The difference between the number of unknowns in 
the system and the number of constraints (or equations) is equal to the number of degrees of freedom 
/ in the system. 


f = (CJ + 1)7T — (CJ + 2) (7 T ~ 1) 

= 2 + cj-tt (3.9) 

This is known as the Gibbs phase rule. It tells us the number of variables / that must be specified in 
order to fix the (intensive) state of the system. 



One generation’s transformation is the next’s status quo. 
In the near future, people may soon think it’s strange that 
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Figure 3.1: Phase diagram for mixtures of methane and ethane. Data taken from a variety of sources. 


3.2 Binary mixtures 

3.2.1 Vapor-liquid equilibrium 

The Gibbs phase rule for a binary mixture simplifies to 

/ = 4-tt (3.10) 

The maximum number of degrees of freedom that a system can have is given when there is only one 
phase present (i.e., n = 1). For a binary mixture, we find that there are at most 3 degrees of freedom. 
This means that we can represent the state of binary mixture using a three dimensional diagram. An 
example of such a diagram is given in Fig. 3.1, which is for mixtures of methane and ethane. 

The key feature of this phase diagram is a solid body in the center of the figure. Within this solid 
body, the system exists as a two-phase mixture, with a coexisting liquid and vapor phases. Above this 
body, the system exists as a single liquid phase; below this body, the system is a single vapor phase. 
The upper surface (marked by red points) that bounds the body is the locus of bubble points (i.e., the 
points at which bubbles begin to appear in a liquid). The lower surface (marked by green points) is 
the locus of dew points (i.e., the points at which droplets begin to appear in a vapor). 
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temperature (K) 


Figure 3.2: Generic pressure-temperature diagram for binary mixtures of methane and ethane: (i) 
pure methane (black line), (ii) 15mol% ethane (red lines), (iii) 50mol% ethane (green lines), (iv) 
70 mol% ethane (blue lines), and (v) pure ethane (violet line). The solid lines and filled symbols 
denote the bubble point curves (saturated liquid), and the dashed lines and open symbols denote the 
dew point curves (saturated vapor). Data taken from RT Ellington et al., Pap. Symp. Thermophys. 
Prop. 1, 180 (1959). 


The points C\ and C 2 are the critical points of pure methane and ethane, respectively. The line 
connecting these two points, which is the intersection of the bubble point and dew point surfaces, is 
the critical locus. This is the set of critical points for the various mixtures of methane and ethane. The 
black curve connecting points A and C\ is the vapor pressure curve of pure methane, and the violet 
curve connecting points B and C 2 is the vapor pressure curve of pure ethane. 

We can represent Fig. 3.1 in a two-dimensional figure by taking various projections. The blue enve- 
lope is a horizontal cross-section of the two-phase body; this is a Txy diagram of the methane-ethane 
mixture taken at constant pressure. The brown envelope is a vertical cross-section of the solid body, 
taken at constant temperature; this is a pxy diagram. 

In Fig. 3.2, we show the pressure-temperature view of the phase diagram for binary mixtures of 
methane and ethane. The point C\ represents the critical point of pure methane, and the point C 2 
represents the critical point of pure ethane. The curve connecting the points A and C\ is the vapor 
pressure curve for pure methane; the curve connecting points B and C 2 is the vapor pressure curve 
for pure ethane. The dotted curve connecting the points C\ and C 2 is the critical locus. The critical 
points of the mixtures, where the coexisting liquid and vapor phases become identical, lie on this 
critical locus. 

For a one component system, the bubble point and the dew point are the same and lie along the vapor 
pressure curve; however, this is not necessarily the case for a mixture. Within envelopes contained 
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mole fraction of toluene mole fraction of toluene 


Figure 3.3: Vapor-liquid equilibria of mixtures of cyclohexane and toluene: (a) Txy diagram at a 
pressure of 1 atm, and (b) pxy diagram at a temperature of 50° C. The filled symbols denote the 
bubble point curves (saturated liquid), and the open symbols denote the dew point curves (saturated 
vapor). Data taken from F Rivenq, Bull. Soc. Chim. Fr. 9, 3034 (1969). 


between the vapor pressure curves of the pure components, a mixture consists of a coexisting vapor 
and liquid phases. The upper part of the envelop (the solid curve with filled symbols) is the bubble 
point curve; the lower part of the envelop (the dashed curve with open symbols) is the dew point 
curve. Different envelopes correspond to different mixture compositions. 

In Fig. 3.3a, we present the Txy diagram for binary mixtures of cyclohexane and toluene at a pressure 
of 1 atm, which is below the critical pressure of both pure species. Point A denotes the boiling 
temperature of pure toluene, and point C is the boiling temperature of pure cyclohexane. Connecting 
these two points are two curves that form the two-phase envelope. The upper curve (with the open 
symbols) is the dew point curve, and the lower curve (with the filled symbols) is the bubble point line. 

Above the two-phase envelope, the system is a vapor, and below the envelope, the system is a liquid. 
Within the envelope, the system separates into a coexisting vapor and liquid phase. The composition 
of the phases is given by the dew point curve and the bubble point curve. For example at point E 
(mole fraction z a \ the system splits into a vapor phase with a composition corresponding to point 
D (mole fraction y a ) and a liquid phase with a composition corresponding to point B (mole fraction 
x a ). The ratio of the total moles of the liquid phase N ^ to the total moles of the vapor phase N ^ is 
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given by the lever rule: 


z a N = x a N ® + y a N 

z a (N (1) + iV (s) ) = x a M l) + y Q iV (s) 

^ (1) Zq - y a 
N(s) x a - z a 

where N = N® + is the total number of moles in the system. 


(3.11) 


In Figure 3.3b, we show the pxy diagram for binary mixtures of cyclohexane and toluene at 50°C, 
which is below the critical temperature of both species. Point A is the boiling pressure of pure toluene, 
and point C is the boiling pressure of pure cyclohexane. Connecting these two points are the bubble 
point (upper) and dew point (lower) curves. Above the bubble point curve, the system is entirely in 
the liquid phase, while below the dew point curve, the system is in the vapor phase. Between these 
two curves, the system separates into a coexisting liquid and vapor phase. The lever rule, described 
previously for the Txy diagram, also applies to the pxy diagram and can be used to determine the 
relative proportions of these phases. 
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mole fraction of methane mole fraction of methane 

Figure 3.4: Phase diagram for mixtures of methane and ethane, (a) Txy diagram with pressure: (i) 
0.101325 MPa (red curves), (ii) 0.6895 MPa (green curves), and (iii) 5.861 MPa (blue curves), (b) 
pxy diagram with temperature: (i) 172.04 K (red curves), (ii) 199.93 K (green curves), and (iii) 280 K 
(blue curves). Data taken from RT Ellington et al., Pap. Symp. Thermophys. Prop. 1, 180 (1959). 


Now let’s consider the variation of the Txy diagram of a binary mixture with pressure. In Figure 3.4a, 
we show the Txy diagram for mixtures of methane and ethane. At moderately low pressures, we have 
the standard Txy diagram, which touches the temperature axis at the boiling temperature of each 
of the pure species. As the pressure increases, the two-phase envelope gradually moves to higher 
temperatures, due to the fact that boiling temperatures increase with pressure. However, when the 
pressure of the system becomes higher than the critical pressure of methane (46.0 bar), the two-phase 
envelope no longer touches the temperature axis at ^methane = 1. As the pressure further increases 
beyond the critical pressure of ethane (48.8 bar), then the two-phase envelope also detaches from the 
temperature axis at x m ethane = 0. 

In the systems that we have examined so far, the bubble point and the dew point of the mixture vary 
monotonically with the composition. This is the case for ideal systems. However, for very non-ideal 
systems, there may be a maximum or a minimum in the bubble and dew point curves. This is the case 
for azeotropic systems. An example of a system that exhibits a low-boiling azeotrope is a mixture of 
n-heptane and ethanol, which is shown in Figure 3.5. For this type of system, both the bubble and 
dew point temperature curves have a local minimum at the same composition. At this composition, 
these two curves meet. This point is known as the azeotrope. At the azeotrope, the composition of the 
coexisting liquid and vapor phases are identical. In this case at the azeotrope, the boiling temperature 
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mole fraction of ethanol mole fraction of ethanol 


Figure 3.5: Phase diagram for mixtures of n-heptane and ethanol, which exhibits a low boiling 
azeotrope: (a) temperature-composition diagram at 30.1°C, and (b) pressure-composition diagram 
at 1 atm. Data taken from JD Raal, RK Code, and DA Best, J. Chem. Eng. Data 17, 21 1 (1972). 


of the liquid is lower than the boiling temperature of either pure components. The corresponding 
bubble and dew point pressure curves have a maximum at the azeotrope (see Fig. 3.5b). 

For a high boiling azeotropic system, the bubble and dew point temperature curves meet at a maxi- 
mum in the Txy diagram. Mixtures of nitric acid and water form exhibit a high boiling azeotrope; 
this system is shown in Fig. 3.6. 


3.2.2 Liquid-liquid equilibria 

When two liquids are mixed together, they do not always form a single, homogenous liquid phase. In 
many cases, two liquid phases are formed, with one phase richer in the first component and the other 
phase richer in the second component. The classic example of a system that exhibits this behavior is 
a mixture of oil and water. 

A typical Txy diagram for a system that demonstrates liquid-liquid phase separation is given in 
Fig. 3.7, which is for mixtures of phenol and water. Outside the curve, the mixture exists as a single, 
homogeneous liquid phase. Inside the curve, the system exists as two separate, coexisting liquid 
phases. Consider the point A , which lies inside the curve. In this case, the system phase separates 
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mole fraction of nitric acid mole fraction of nitric acid 


Figure 3.6: Phase diagram for mixtures of nitric acid and water, which has a high boiling azeotrope: 
(a) temperature-composition diagram at 1 atm, and (b) pressure-composition diagram at 25°C. Data 
taken from MA Yakimov and VYa Mishin, Radiokhimiya 6, 543 (1964). 


into two liquid phases, one with composition x[ and the other with composition x'[. The relative 
amounts of the two phases can be computed using the lever rule, as described previously. 

The maximum of the liquid-liquid phase envelope and is known as the critical point of the mixture. 
Above the critical temperature (i.e., the temperature at the critical point), the system exists as a single 
liquid phase. Below the critical temperature, the system can split into two coexisting liquid phases, 
depending on the overall composition. 

The basic reason why liquid-liquid phase separation occurs is that the attractive interactions between 
different molecules are weaker than the attractive interactions beween similar molecules. As a result, 
similar molecules prefer to be near to each other and than to dissimilar molecules. The consequence 
of this at a macroscopic level is the formation of two liquid phases. 


3.2.3 Vapor-liquid-liquid equilibria 

In Figure 3.8, we combine the high temperature and low temperature Txy phase diagrams for a system 
with a low boiling azeotrope. At high temperatures, the system is a vapor. At lower temperatures, 
it condenses to become a liquid. Below the critical temperature, liquid can separate to form two 


31 


Download free eBooks at bookboon.com 


Chemical Thermodynamics 


Multicomponent systems 



Figure 3.7: Phase diagram of mixtures of phenol and water. Data taken from AN Campbell and A JR 
Campbell, J. Am. Chem. Soc. 59, 2481 (1937). 


coexisting liquid phases. 

As the pressure of the system decreases, the boiling temperature decreases, in general. Therefore, we 
expect the vapor-liquid coexistence envelope to drop to lower temperatures as the pressures decrease. 
Changes in pressure, however, do not have a strong influence on the phase behavior of liquids. As a 
result, we do not expect the liquid-liquid phase envelope to change much with pressure. Consequently, 
we expect that at a low enough pressure, the vapor-liquid coexistence curve will intersect the liquid- 
liquid coexistence curve. When this occurs, we can have vapor-liquid-liquid equilibria, which is 
shown in Figure 3.9. 
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Figure 3.8: Phase diagram for system with a low boiling azeotrope and two liquid phases: (a) 
temperature-composition diagram, and (b) pressure-composition diagram. 




Figure 3.9: Phase diagram for system with a heterogeneous azeotrope: (a) temperature-composition 
diagram, and (b) pressure-composition diagram. 
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1 



Figure 3.10: Triangular phase diagram representing a ternary mixture. The filled circle represents a 
40%, 40%, 20% mixture of components 1, 2, and 3, respectively. 


3.3 Ternary mixtures 

In this section, we consider the phase behavior of systems with three components. The Gibbs phase 
rule for a ternary mixture is 

/ = 5-tt (3.12) 

A ternary system can have at most four degrees of freedom, which occurs when there is only one 
phase present (i.e., n = 1). At a given temperature and pressure, only the composition of a ternary 
system needs to specified to fix its state. 

In this case, the phase diagram is typically represented on an equilateral triangle, such as the one 
shown in Figure 3.10. The comers of the triangles represent the pure components. The upper comer 
is pure 1; the lower left comer is pure 2; and the lower right comer is pure 3. For a general point 
within the triangle, the mole fraction of a species is equal to the ratio of the distance of the point from 
the corresponding comer to the height of the triangle. Therefore, the sides of the triangle represent 
binary mixtures. 

Some typical phase behavior that can be exhibited by ternary mixtures is shown in Fig. 3.11. Let 
us consider a situation where binary mixtures of component 1 and component 2 are only partially 
miscible, where two coexisting liquid phases may be formed: one rich in 1 and the other rich in 
2. This is represented by the base of the ternary phase diagram shown in Fig. 3.11a. In addition, 
let us assume that components 1 and 3 are completely miscible and components 2 and 3 are also 
completely miscible. For this case, one might expect that if enough of component 3 is added to the 
system, then components 1 and 2 can be made to mix with each other, due to their mutual solubility 
with component 3. This is type I phase behavior. 
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Figure 3.1 1: Generic phase behavior of ternary mixtures: (a) type I, (b) type II, and (c) type III. 


A type II phase diagram, shown in Fig. 3.11b, corresponds to a situation where components 1 and 
3 are completely miscible, but both components 1 and 2 and components 2 and 3 are only partially 
miscible. 


Finally, type III phase behavior is shown in Fig. 3.11c. In this case, the various binary mixtures of 
the three components are each only partially miscible. The shaded triangle in the center of the phase 
diagram is a region where three phases are in coexistence with each other. Systems with a composition 
which lies within this shaded triangle will split into three separate phases; the composition of each 
of these phases corresponds to one of the comers of the triangle. The composition of the individual 
phases will not vary with the system’s location within the triangle (i.e., its overall composition); 
however, the relative amounts of each of the phases will. 
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4 The ideal solution model 


In many situations, we need to predict the properties of a mixture, given that we already know the 
properties of the pure species. To do this requires a model that can describe how various components 
mix. In mathematical terms, this means that we need to relate the Gibbs free energy of a mixture to 
the Gibbs free energy of the various pure components. One of the simplest models that achieves this 
is the ideal solution model. In this lecture, we present the ideal solution model. Then we apply this 
model to describe vapor-liquid equilibria, and as a result, derive Raoult’s law. 


4.1 Definition of the ideal solution model 


For an ideal solution, the Gibbs free energy G defined as: 

G(T, p,N 1 ,N 2 ,...) = J2 N afi ° a (T, p) + RTj2 N a In (4.1) 

a a 

where T is the absolute temperature of the system, p is the pressure of the system, R is the ideal 
gas constant, N a is the number of moles of component a, x a is the mole fraction of component a, 
and /i° (T, p) is the molar Gibbs free energy of pure component a (recall that for pure systems, the 
molar Gibbs free energy is equal to the chemical potential). The first term on the right-hand side of 
Eq. (4.1) represents the Gibbs free energy of the system if its components were unmixed. The second 
term represents the contribution due to the entropy of mixing. 

In all calculations involving the ideal solution model, we assume that we know the molar Gibbs free 
energy of each of the pure species as a function of temperature and pressure. Mathematically, this 
means that know the form of the functions /i° (T,p). Physically, this means that we know everything 
about the thermodynamics of the pure species. 


Once we know the Gibbs free energy of a system as a function of temperature, pressure, and compo- 
sition, we know everything about its thermodynamics. For example, the total volume of the system 
Y can be derived from the Gibbs free energy: 


F=r f) 


T,N a 




T,N a 




(4.2) 


where V° is the molar volume of pure species a , and we have used the relation V° = (dp^/dp). 
The total volume is equal to the sum of the volumes of the pure components. Therefore, there is no 
change of volume on mixing for an ideal solution. 
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The total entropy S of the system is given by 



= 2 N a S° a (T,p ) - R N a lnx Q 

a a 


(4.3) 


where 5° is the molar entropy of pure component a, and we have used the relation S° = — (dju^ / dT ) . 
Unlike the case for volume, the total entropy of the mixture is not the same as the sum of the entropies 
of the unmixed systems. In fact, due to the second term in Eq. (4.3), the entropy increases upon mix- 
ing. 


The total enthalpy H of the system is given by 


h = g + ts 


= ]T iV Q /i° (T, p) + RT ]T N a In + T I £ N a S° a (T, p) 

a a \ a 

= '£N a H° a (T,p) 

a 


ln 


X n 


(4.4) 


where is the molar enthalpy of pure a , and we have used the relation = /i° + TS°. In this 
case, we see that the total enthalpy of the system is the same as the sum of the enthalpies of the pure 
components. No heat is absorbed or released upon mixing an ideal solution. 


The chemical potential of species a is given by 


[lex, (-^ Pi X2 5 X3 , . . . ) 


\ dN ° ) T ,p, Na r 


= Pa(T, p) + RT In x a 


(4.5) 


In an ideal solution, we see that the chemical potential of a species depends on its mole fraction and 
not directly on the composition of the other components in the system. Also, we see that mixing 
causes the chemical potential of each component to decrease. 


4.2 Derivation of Raoult’s law 

Now we will use the ideal solution model to develop a mathematical description of vapor-liquid 
equilibrium in a multicomponent solution. We will make the assumption that we have a system that 
is separated into a coexisting vapor and liquid phase. The vapor phase will be assumed to behave like 
an ideal gas, while the liquid phase will be assumed to behave as an ideal solution. 
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The basic condition for equilibrium between phases is 

V V a (T,P,y2,y3,---) = Pa(T,p,X2,X 3 ,...) (4.6) 

where /i^ is the chemical potential of component a in the vapor phase, /J a is the chemical potential 
of component a in the liquid phase, y a is the mole fraction of component a in the vapor phase, and 
x a is the mole fraction of component a in the liquid phase. 

Since the liquid phase behaves as an ideal mixture, we have 


Pa (T, P, X 2 , X 3 ,...)= /x°’' (T, p) + RT In x a 

= ^\T,p™v)+[ P d p' d -^^ + RT\nx a 

Jp ™ p op 

= p^iry^) + r dp'V°\T,p') + RT In x a (4.7) 

where we have used the relation {dfi°a /dp) = Va l , and V^’ 1 is the molar volume of pure a in the 
liquid phase. If we assume that the volume of the liquid is nearly independent of pressure, then 
Eq. (4.7) becomes: 


p l a(T,p, X 2 ,x 3 , ...)= p°d\Ty*v) + V°’ l (T,p)(p - P™V) + RT In x a (4.8) 

Since the vapor phase can be considered an ideal gas, we have 

P v a (T,p, 2 / 2 , 2 / 3 , • • • ) = Pa V ( T iP ) + RTlny a 

= p7(.T,pD+ I* dpl 0i ^^ f) + RT\ny a 

^Pa aP 

= P°X(T,pD + r dp'V°’ v (T,p') + RT In y a (4.9) 

•/ Pa aP 

where Va’ v is the molar volume of pure a in the vapor phase. For an ideal gas, we have Va’ v (T, p) = 
RT Ip, and, therefore, 

r p RT 

pl(T,p,y 2 ,y 3 ,...) = ^(T,p™v)+ / dp'— + RT\ny a 

Jp™ p P 

= M°’ v T,pD + RT ln4p + RTlny a (4.10) 

Pa 


Inserting Eqs. (4.8) and (4.10) into Eq. (4.6) yields 


P°a V (T, P ™ p ) + RT In 4p + RT In y a = p° a \T,p^) + V°’ 1 (T,p)(p - C P ) + RT\nx 0 


vap 

Pa 


(4.11) 


We note that at the vapor pressure of a pure system, the chemical potentials of the liquid and vapor 
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phases are equal. That is, /ia V (T,p™ p ) = /i^ 5l (T, j;^ ap ). Substituting this relation into Eq. (4.11) 
gives 


P 


RTln-^p+ RT In y a = V°’\T, p) {p - p™ p ) + RT In 

Pa 


VaP = x a pl ap exp 


RT 


(4.12) 


The exponential term on the right-hand side of Eq. (4.12) is known as the Poynting factor. For most 
systems at low to moderate pressures, the Poynting factor is almost equal to one. In the case of water 
at 25°C: 


Va(P~P V * P ) _ (10“ 6 m 3 /g)(18.02 g/mol)(101325 - 3170) Pa 
RT ~ (8.314 J mol -1 K -1 )(298.15 K) 

» 7.14 x 10 -4 


(4.13) 


This leads to a Poynting factor of 1.000714, which is essentially equal to one. If we assume that the 
Poynting factor is close to one, we have Raoult’s law: 


VaP = X a p Y a V i T ) 


(4.14) 
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5 Partial molar properties 

In Chapter 4, we examined the properties of ideal solutions. Many properties of an ideal solution do 
not change on mixing. For example, the volume of a mixture is equal to the sum of the volume of 
the original unmixed solutions. In this situation, it is straightforward to assign how much volume is 
occupied by each component in the system — it is simply the volume occupied by components in 
their unmixed state. 

For a general system, however, the volume, as well as other properties, is not additive. That is, the 
volume of a mixture is not equal to the sum of the volumes of the individual pure components. In 
this situation, it is not clear how to assign how much volume is occupied of each species. One logical 
manner to do this is through the use of partial molar properties. 

In this Chapter, we define partial molar properties and describe their application. We then discuss 
their relationship with the change of properties of a system on mixing. Finally, we examine the 
graphical representation of partial molar properties for binary mixtures. 

5.1 Definition 

In general for any extensive property X of a system, we define a partial molar property of component 


a as 



(5.1) 


where N a is the number of moles of species a , T is the temperature, and p is the pressure of the 
system. Physically, the partial molar quantity X a corresponds to the change of the property X with 
the addition of a small amount of component a, while holding constant the temperature, pressure, and 
number of moles of all other species. Note that partial molar properties are intensive and, therefore, 
do not depend on the system size. 

Examples of partial molar properties include the partial molar enthalpy H a , which is defined as 



(5.2) 


and the partial molar volume V a , which is defined as 



(5.3) 


Note also that from the relation (see Eq. (3.4)) 



(5.4) 
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we see that the chemical potential is equal to the partial molar Gibbs free energy. 


5.2 Relationship between total properties and partial molar properties 

Any extensive property of a system can be written in terms of its partial molar quantities. Recall that 
an extensive property is a property that scales proportionally with the size of the system. If the system 
doubles in size, then value of the extensive property should double. If the size of the system increases 
by a factor t, then the value of the extensive property should increase by a factor t. For example, 
taking X to be the volume V, then we expect the total volume of the system V to increase by a factor 
t if the number of total moles in the system are increased by a factor £, holding the composition of the 
system fixed. This feature can be expressed mathematically as 

£X(T,p,TV l5 TV 2 , • • • ) = X(T,p,tN u tN 2 , • • • ) (5.5) 


Taking the derivative of both sides of the Eq. (5.5) with respect to t, while keeping all other variables 
constant, yields 



= J]X Q iV Q 

a 


X(T,p,N 1 ,N 2 ,...) = J2 N «X« (5-6) 

a 

The value of any extensive property of a system is equal to the sum of the partial molar properties 
of each component multiplied by the amount of each component in the system. Therefore, we can 
“divide” the property of a mixture , such as the volume or enthalpy, between its individual components 
according to their partial molar properties. 

One important example of this, which we will utilize later, is the Gibbs free energy G. In this case, the 
total Gibbs free energy G can be expressed in terms of the chemical potentials /i a of each component 
(partial molar Gibbs free energy) 


G(T,p,N ll N 2 ,...) = Y,NaVa 

a 


(5.7) 


By dividing both sides of Eq. (5.6) by TV, the total number of moles in the system, we find that molar 
properties are similarly related to partial molar properties 

X = J2 X »Xa (5.8) 

a 


where x a is the mole fraction of component a. 
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The derivative of a molar property with respect to x a can also be written in terms of partial molar 
properties. To demonstrate this, let us consider that an extensive property X is a function of T, p, 
mole numbers N a . We can express its differential as 


dX = 



(5.9) 


Now, we take the derivative of X with respect to x a , holding T, p, all other mole fractions (with the 
exception of x\\ and total number of moles N constant: 


dX 




r)rf I / J ^ \ Fip 

T,p,x a ,^ a l ,N a , V “ / T,p,x a ,^ a l ,N 


dN n 


N 


dX 

dx a 

dX 

dx a 


= N(X a - X 1 


T ’P’ x c,'jtoc,l 


= x a -x 1 


(5.10) 


i 

where we have used the relations N a = Nx a and x\ = 1 — x »- 
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5.3 Properties changes on mixing 

Often when pure components are mixed together to form a solution, the properties of the overall 
system change. For example, if we mix a quantity of pure liquid ethanol with pure liquid water, the 
volume of the final mixture has a slightly smaller volume than the sum of the volumes of the original 
pure liquid. This system has a negative volume change on mixing. If a solution of concentrated 
sulfuric acid is mixed with pure water, then a large amount of heat is released because the enthalpy of 
the mixture (at equilibrium) is lower than the sum of the enthalpy of the unmixed components. This 
system has a negative enthalpy of mixing. 

To characterize this change, the property change on mixing AX mix is defined as the difference be- 
tween the property of the mixture and the properties of the pure components: 

AX mix = X N aX° 

a 

= Y J N a {X a -X° a ) (5.11) 

a 

where X° is the molar property of pure a. If the property increases on mixing, then AJ mix is 
positive; if the property decreases on mixing, then AX mix is negative. 

Similarly, the molar change of a property on mixing AX mix is given by: 

A X mix = J2*a(Xa-X°) (5.12) 


5.4 Graphical representation for binary systems 

The relations derived in this section provide us a means to easily extract the partial molar properties 
of a system from a graph of the corresponding molar property with repect to composition. A generic 
plot of the variation of a molar property of a binary mixture with composition is given by the solid 
line in Fig. 5.1. When X 2 = 0, the system consists only of component 1, and the value of the molar 
property should be equal to the molar property of pure component 1 X%. Likewise, when X 2 = 1, 
the system consists only of component 2, and the value of X should be equal to X$. The variation of 
a molar property of an unmixed system composed of pure 1 and pure 2 is given by the thick dashed 
line in Fig. 5.1. 

The distance between the thick dashed line and the solid line represents the difference between a 
molar property of the mixture and that of the unmixed system. This difference is equal to AX mix . 

Now consider a line that is tangent to X (i.e., the solid line) at the composition X 2 = x\. This is 
denoted by the thin dashed line in Fig. 5. 1 . The intercept of the tangent line at X 2 = 0 and X 2 = 1 are 
equal to X\ and X 2 , respectively; that is, the intercepts of the tangent line are equal to the partial molar 
properties of the system at a given composition. To demonstrate this, let’s determine the equation for 
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X 2 

Figure 5.1: Variation of the molar property X (e.g., molar enthalpy or molar volume) of a binary 
mixture as a function of the mole fraction of component 2 (solid line). The dashed line is the tangent 
line to X at the composition X 2 = x° 2 . 


the tangent line y[x 2 ) at x^ = x\. 

At x 2 = x 2 , the tangent must have the value X{x 2 )\ 

y(x 2 ) = 

= ^X!(^)+^X 2 (^) 

= XiixZ) + x° 2 (X 2 (x° 2 ) - X^xD), (5.13) 

where X a (x 2 ) is the value of the partial molar property of component o at x 2 = x 2 , and we have 
used Eq. (5.8). Note that the partial molar properties X\ and A’_> are usually functions of composition 
(i.e., x 2 ). The slope of the tangent line can be obtained from Eq. (5.10) 

{^) T , r = Ux%) ~ Uxl) - <5 ' 14) 

Therefore, the slope of the tangent at composition X 2 = x° 2 is equal to X 2 ^X 2 ) — X\{x 2 ). 

Knowledge of both the value of the tangent line at a point and its slope allows us to write the equation 
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for the tangent line: 

y(x 2) = y(x 2) + [x 2 (x%) - X 1 {x 0 2 )\{x2 - x° 2 ) 

= X^xl) + x° 2 (X 2 (x° 2 ) - X^xD) + [X 2 (x° 2 ) - Xi(x° 2 )](x 2 - x ° 2 ) 

= X 1 (x 0 2 ) + [X 2 (x° 2 )-X 1 (x° 2 )}x 2 (5.15) 

From Eq. (5.15), we see that when X 2 = 0, the tangent line has a value equal toXi^), and when 
X 2 — 1? it has a value equal to This is shown graphically by the thin dashed line in Fig. 5.1. 

For an ideal solution, the molar volume V of a mixture satisfies 


V = x 1 y l ° + x 2 y 2 ° (5.16) 

where V x ° and y 2 ° is the molar volume of pure species 1 and 2, respectively. Therefore, Al/ mix = 0 
for an ideal solution. Similarly, the molar enthalpy H is given by 

H = x x H{ + x 2 H ° 2 (5.17) 

where H\ and is the molar volume of pure species 1 and 2, respectively. Therefore, AiT mix = 0 
for an ideal solution. 
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Equations (5.16) and (5.17) correspond to a straight line connecting the molar properties of the pure 
components (e.g., the thick dashed line in Fig. 5.1). In this case, we see that the partial molar volumes 
and partial molar enthalpies of each species are equal to the respective molar quantities in the pure 
state. We expect ideal behavior when the fluids that are mixed consist of similar molecules. 

Most solutions do not exhibit ideal behavior, and the actual curve corresponding to the variation of the 
molar volume or enthalpy of the mixture deviates from a straight line (e.g., the solid line in Fig. 5.1). 
When the curve for the molar volume lies above the ideal mixture line, the system expands upon 
mixing; when the curve lies below the line, the system contracts. In the case of the molar enthalpy, a 
curve that lies above the ideal mixture line corresponds to the system that absorbs heat (e.g., mixing 
lead bromide and water); a curve that lies below the line corresponds to the system releasing heat 
(e.g., mixing sulfuric acid and water). This non-ideal mixing in the case of the molar enthalpy is 
the principle used in cold packs and heat packs. We will develop mathematical models to describe 
non-ideal mixtures. We use partial molar properties in more detail later. 
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6 Nonideal solutions 

In Chapter 4, we developed the ideal solution model, which enables the estimation of the properties 
of mixtures from knowledge of the thermodynamic behavior of the pure species. While the ideal 
solution model does provide accurate predictions for mixtures of relatively similar substances, many 
systems do exhibit substantial deviations from the ideal solution model. 

In this chapter, we present methods for mathematically describing the properties of non-ideal solu- 
tions — mixtures that deviate from the ideal solution model. 


6.1 Deviations from Raoult’s law and the activity coefficient 

As discussed in Chapter 4, the vapor-liquid equilibrium between an ideal gas and an ideal liquid 
solution is governed by Raoult’s law: 


y a p = x a p™ p (T), 

where y a is the mole fraction of component a in the vapor phase, x a is the mole fraction in the liquid 
phase, p is the pressure of the system, and p™ p (T) is the vapor pressure of pure component a at the 
system temperature T. 

Raoult’s law states that the partial pressure of a species in the vapor phase is proportional to the mole 
fraction of the species in the liquid phase. The constant of proportionality is the vapor pressure of the 
pure species. Mathematically, this expressed as 

Pa = X a p™ P (T). (6.1) 

This equation is plotted as the solid line in Fig. 6.1. 

For most systems, Raoult’s law is not obeyed. If we actually measure the partial pressure of a given 
component above a mixture as a function of its mole fraction, it will most often deviate from a straight 
line. For example, see the dashed and dotted lines in Fig. 6.1. In order to characterize the deviations 
from Raoult’s law, we define the activity coefficient for component a , which is denoted by as: 

7a = zrfVr (6 ' 2) 

Note that in general, the activity coefficient depends on the temperature, pressure, and composition 
of the system. If the system obeys Raoult’s law, = 1. When > 1, the system is said to 
exhibit positive deviation from Raoult’s law, and when < 1, the system is said to exhibit negative 
deviation from Raoult’s law. 

Although Raoult’s law is not a good approximation for many systems, there are two points at which 
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^ct 

Figure 6.1: Variation of the partial pressure of component a in the vapor phase with its mole fraction 
in the coexisting liquid phase. The solid line is for a component that obeys Raoult’s law, the dashed 
line is for a component that shows positive deviations from Raoult’s law, and the dotted line is for a 
component that shows negative deviations from Raoult’s law. 


it is exact: (i) at x a = 0, where the partial pressure of a must equal zero, and (ii) at x a = 1, where 
the partial pressure of a must be equal to p™ p (T). From condition (ii), combined with Eq. (6.2), the 
activity coefficients must obey: — ► 1 as x a — > 1. 


6.2 Modified Raoult’s law 

As mentioned previously, many systems show significant deviations from Raoult’s law. There are 
even some solutions, such as azeotropic mixtures, that are not even qualitatively described by Raoult’s 
law. By rearranging the definition of the activity coefficient (see Eq. (6.2)), Raoult’s law can be 
generalized to apply to a wider class of mixtures: 

VaP = XalaP V a V { T )- ( 6 - 3 ) 

This equation is known as the modified Raoult’s law. 

Raoult’s law was derived under the assumptions that the vapor phase behaved as an ideal gas and 
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that the liquid phase behaved as an ideal solution. In this section, we will relax the assumption that 
the liquid behaves as an ideal solution to derive the modified Raoult’s law, given in Eq. (6.3). If we 
re-examine the derivation of Raoult’s law, we find that we can arrive at the modified Raoult’s law if 
the chemical potential of the liquid phase can be written as: 

Z-4 x 2 , x 3 , . . . ) = (T, p) + RT In x a ^ a (6.4) 

where n°d l (T, p) is the chemical potential of pure component a. Note that we can divide the chemical 
potential into an ideal contribution /i 1 ^, given by the ideal solution model, and a non-ideal contribution 
(or excess) /i^ x , which is directly related to the activity coefficient: 

/i a (T,p,x 2 ,x 3 ,...) = /i^(T,p,x 2 ,x 3 ,...) +^ x (T,p,x 2 ,x 3 ,...), (6.5) 


where 


Va(T,P, x 2 ,x 3 ,...) = n°J(T,p) + RTlnx a , (6.6) 

and 

fj, e *(T,p,x 2 ,x 3 ,...) = RTln'y a (T,p,x 2 ,x 3 ,...). (6.7) 

In general, an excess property X ex is defined as the deviation of a property from its ideal solution 
value: 

X ex (T,p, X 2 , X 3 ,...)= X{T, p,x 2 ,x 3 , . . . ) - X id (T,p, X 2 , x 3 ,...) (6.8) 

Substituting the expression for the chemical potential given in Eq. (6.4) into Eq. (5.7), we find that 
the Gibbs free energy of a non-ideal solution can be written as: 

G — ^ X a fi a 

a 

= ^2 N a (Xa + RT In X a 7a) 

a 

= J2 N M + RTlnx a ) + RTJ2 Na In 7a 

a a 

G = G id + G ex (6.9) 

where G ld is the Gibbs free energy of an ideal solution, defined as 

G id = ^ /V qM ° + RT J2 Na In a; Q (6.10) 

a a 

and G ex is the excess Gibbs free energy, which is given by 

G ex = RT N a In 7 q (6.11) 
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While the theoretical treatment we have presented for non-ideal solutions are all formally exact, they 
are of no practical use until we have some kind of model for the variation of the activity coefficients 
with the system composition, temperature, and pressure. In the next section, we present various 
models that can be employed to describe the activity coefficients. 
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6.3 Empirical activity coefficient models 

In order to perform quantitative thermodynamic calculations using the Gibbs free energy for a non- 
ideal solution (see Eqs. (6.9)-(6.1 1)), we need explicit expressions for the activity coefficients. A few 
empirical expressions that are typically employed are: 

• two-suffix Margules 

In 71 
In 72 

G cx 

RT 

• three-suffix Margules 

ln7i = [(M + 3 B) — 4Bx 2 }x\ 

In 72 = [(M - W) + 4Bxi]xj 

£r GX 

— = [A + B(xi - x 2 )]x 1 x 2 

• four-suffix Margules 

In 71 = [( A + 3 B + 5 C) — 4 (B + 4C)x2 + VlCx^x\ 

In 72 = [( A — 3 B + 5 C) + 4 (B — AC)x\ + 12 Cx\\x\ 

G ex 

— = [A + B(xi - x 2 ) + C(x i - x 2 ) 2 ]xix 2 


= Axl 
= Axf 

~ Ax\X2 


• van Laar model 


lnyi — A 
In 72 — B 



2 

Bx 2 ) 

Bx2_y 2 

Ax i ) 


G ex _ ABx ix 2 
RT Ax i + Bx 2 


The parameters A, B , and C in the above models are, in general, functions of temperature and 
pressure, but are independent of composition. These parameters are typically obtained by fitting 
experimental data. Given the parameters of these activity coefficient models, we can predict all the 
thermodynamic properties of the system. 
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6.4 The Gibbs-Duhem equation 

Not all the chemical potentials (and therefore, the activity coefficients) in a mixture are independent 
of each other. They are all related to one another through the Gibbs-Duhem equation. To derive this 
equation, we start with the fundamental equation of thermodynamics for the Gibbs free energy, which 
can be written as 


dG = —SdT + Ydp + Y VadN a (6.12) 

a 

By substituting Eq. (5.7) into Eq. (6.12), we find 

d NaUc^j 

Y + Y ^ dN « 

a a 

Y N a d l l a 

a 

This relation is known as the Gibbs-Duhem equation. Dividing both sides of the equation by the total 
number of moles in the system, we find that the Gibbs-Duhem equation can also be written as 

x a d[i a = —SdT TV dp (6.14) 

a 

The Gibbs-Duhem equation provides a relation between the chemical potentials of each of the chem- 
ical species in a given phase. At constant temperature and pressure, this relation simplifies to 

Y Xgdpa = 0 (6.15) 

a 


—SdT + Ydp + p a dN a 

a 

-SdT + Ydp + Y HadN a 

a 

-SdT + Ydp (6.13) 


The chemical potential can be written in terms of the activity coefficient (see Eq (6.4)). Subsituting 
this expression into the Gibbs-Duhem equation, at constant temperature and pressure, we find 


Y X ad(/J,° a + RT In x a j a ) 

a 

RT E x a d\n XoTia 

a 


RT x a (d lnx a + dlny a ) 

a 

Y Xa (— + din 7«) 

a \ Xa ' 

Yjydx a + Xadln^a) 

a 

Y x »d In 7a 

a 


o 

o 

0 

0 

0 

0 


(6.16) 
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where in the last line we have used the relationship d x = 0 . This equation provides a rela- 
tionship between the activity coefficients of all of the components in a system. Therefore, not all 
the activity coefficients in a system are independent. If we know the activity coefficients for all but 
one of the components in a mixture, we can determine the unknown activity coefficient by using the 
Gibbs-Duhem equation. We demonstrate this for a binary system. 

Suppose we know the activity coefficient of component 1 71 as a function of composition. Then, we 
can determine the activity coefficient of component 2 72: 


xidlnji + £2<iln72 = 0 

din 71 , din 72 
X \ — b X2 — = 0 


dx 2 


/ X 2 

dx 2 


x\ din 71 
X2 dx 2 


f 

In 72(^2) - In 72(2:2 = 1) = - / 

f 

In 72 {x 2 ) = - J 


dx 2 
d In 72 
dx 2 
d In 72 

dx 2 J 1 ""~*X2 dx 2 

1 — X2 d In 71 


-L 


X2 xidlnji 

dx 2 - 


dx 2 - 


dx 2 


X2 dx 2 
1 — X2 d In 71 
X 2 dx2 


( 6 . 17 ) 


where we have used the fact that 72 — ► 1 as X2 — ► 1 . Given the 71(^2), Eq. ( 6 . 17 ) can be used to 
determine 72(^2)* 

For example, if we are given In 71(2:2) = Ax\ (the two-suffix Margules model), then 

/ *^2 1 HTo d 

dx 2— — — (Ax 2) 


-1; 

= -2 A 
= -2 A 
= - 2 A 


dX 2 


X2 dX2 
1 - X2 


X 2 


( 2 Ax 2 ) 


/ X2 

dx 2 ( 1 - x 2 ) 


XZ 


x 2 


X2 


J 1 


(X 2 ~ 1 ) 


(*2 - !) 


= A( 1 - x 2 f 


( 6 . 18 ) 


6.5 Azeotropic systems 

One of the shortcomings of Raoult’s law is that it cannot describe the behavior of azeotropic systems. 
The modified Raoult’s law, however, can describe azeotropic systems as we will demonstrate below. 
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For a binary mixture, the modified Raoult’s law can be written as 

vap 

yip = 71Z1P! p 

V2P = ^2X2pT P (6.19) 


At the azeotrope, the vapor composition is the same as the liquid composition: 

:C eo = -< zeo (6.20) 

At this point, we can cancel the compositions from both sides of Eq. (6.19), which gives 


P = 71 P? P (T) 

P = l2 p? P (T) (6.21) 

These two equations can be combined into a single equation that can be used to solve for the azeotropic 
composition x| zeo : 


7 i pT p (t) = 12 pT p (t) 

71 = P? P (T ) 

72 


In — = In 

72 


P? P (T) 
PT{T) 
pT p (t) 


( 6 . 22 ) 
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Given an activity coefficient model, we can use Eq. (6.22) to determine the azeotropic composition. 
For example, if we use a two-suffix Margules model, where 

In 7 i = Ax 2 
In 72 = Ax\ 


we have 

In — = — Ax\ 

72 

= A(X2 + Xi)(x2 ~ Xi) 

= A(x 2 ~ Xi) 

= ,4(2x 2 - 1) (6.23) 


By substituting Eq. (6.23) into Eq. (6.22), we find 


,4(2xf eo 


— 1 ) = In 


pT(t} 

pT p (t) 


azeo 

x 2 


1 

2 



pT p ( t ). 


(6.24) 


Therefore, if we are given the parameter A, we can determine the azeotropic composition. 


On the other hand, if we are given the azeotropic composition, then we can fit the parameter A. 
Rearranging Eq. (6.24) yields 


A = 


1 

2xf eo - 1 


In 


PT(T) 

P? P (T) 


(6.25) 


We can also determine what values of the parameter A will lead to an azeotropic system. By noting 
the restriction 0 < x| zeo < 1, a system will have an azeotrope if the A parameter satisfies: 


-In 


pT(T) 

pT p (t) 


< 


A 


< In 


pT\t) 

pT p (t) 


(6.26) 


where we have made the assumption that p™ p (T) > p\ ap (T). 
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7 Stability 

In this Chapter, we consider the mathematical and geometric properties of the free energy of a system 
which govern its thermodynamic stability. In regions where the free energy satisfies the conditions 
for stability, the system remains in a single phase; however, if these conditions are violated, then the 
system can phase separate into two or more coexisting phases. We then discuss how the conditions 
for stability on the free energy leads to conditions on the value of the parameters in activity coefficient 
models. Finally, we derive the mathematical equations for the coexistence between two liquid phases. 


7.1 Introduction 


A system at constant temperature, pressure, and mole numbers tries to decrease its Gibbs free energy. 
At equilibrium, the Gibbs free energy of the system is minimized. 

In Fig. 7.1, the solid line represents the variation of the molar Gibbs free energy of a mixture as a 
function of composition. The thick dashed line represents the variation of the molar Gibbs free energy 
of a completely unmixed system, with one phase consisting of pure component 1 and the other phase 
consisting of pure component 2. The distance between the solid line and the dashed line denotes the 
Gibbs free energy change on mixing. Because the solid line always lies beneath the dashed line, the 
system possesses a lower Gibbs free energy when mixed than when separated and, therefore, prefers 
to be mixed. If the opposite were true, then the system would prefer to be phase separated. 


The solid line in Fig. 7.2 represents the molar Gibbs free energy of mixture that becomes unstable. 
The thick dashed line, as in Fig. 7.2, represents the molar Gibbs free energy of a completely unmixed 
system. The solid line always lies beneath the dashed line, and so the mixture will not phase separate 
to form two completely demixed phases. However, consider the line joining the points A and B , 
which lies beneath the solid line. The line AB represents the molar Gibbs free energy of a system 
that is separated into two phases, one with composition x 2 and the other with composition x 2 . This 
partially demixed system has a lower Gibbs free energy than the fully mixed system (solid line), and, 
therefore, a mixture which is at a composition between x 2 and x 2 will spontaneously demix to form 
these two, coexisting phases. Another point to note, is that because the tangent line at x 2 and x 2 
coincide, the chemical potentials in each of the two corresponding phases are identical, which is the 
criterion for phase equilibrium. 


In general, a mixture is unstable when 


d 2 G 

dx 2 2 


< 0 


(7.1) 


Graphically, this means that the Gibbs free energy is concave down with respect to X 2 (e.g., see the 
solid line in Fig. 7.2 between x 2 and x 2 ). In this situation, any line that connects two points on the 
curve will lie below the curve, and thus the mixture will want to phase separate. 
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Figure 7.1: The molar Gibbs free energy G of a stable binary mixture as a function of the mole 
fraction of component 2 (solid line). The dashed line is the tangent line to G at the composition 

x 2 = x° 2 . 


An ideal solution is never unstable and will always remain as a single phase mixture. To prove this, 
we need to demonstrate that d 2 G ld /dx \ is always positive. For a binary mixture, the molar Gibbs 
free energy of an ideal solution is given by 

G ld — xi/i\ + X2H2 + RTx 1 lnxi + RTx 2 lnx 2 

= (1 - x 2 )/i 1 + x 2 /i 2 + RT( 1 - X2) ln(l - X2) + RTx 2 In x 2 (7.2) 


Taking the derivative of this, we find 

d G* d 

— — = H 2 — /ii — RT ln(l — X 2 ) + RT In £2 

d 2 G id _ RT RT 
dx 2 1 — X2 X2 

RT 

X 2 (l - X 2 ) 


(7.3) 


For all physical values of x 2 (i.e., 0 < .>'■> < 1), we see that d 2 G ld /dx 2 > 0. Therefore, the ideal 
solution model is always stable. 
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X 2 

Figure 7.2: Variation of the molar Gibbs free energy G of a binary mixture as a function of the mole 
fraction of component 2 (solid line). The long-dashed line is the Gibbs free energy for a system 
of two unmixed pure compounds. The short-dashed line is the tangent line to G at the coexistence 
compositions X 2 = x 2 and X 2 = x 2 . 


To model a mixture that phase separates into two coexisting liquid phase, we need to add non-ideal 
terms (activity coefficients) to the ideal solution model. As an example of this, we examine the 
stability of the two-suffix Margules model, which has a molar Gibbs free energy of 

G = x\ Hi + X2H2 + RTx 1 lnxiyi + RTx 2 In £272 

= (1 — X2)/ii + X21R2 + RT( 1 — X 2 ) ln(l — X 2 ) + RTx 2 In £2 + RTA (1 — X 2 )x 2 (7.4) 


The derivatives of the molar Gibbs free energy are: 


dG 

— — = — fi T jji 2 — RT Yn( 1 — X2) T RT In X2 T RTA{\ — 2 ^ 2 ) 


d 2 G 

dx 2 2 


RT RT 

1 2 RTA 

1 - x 2 x 2 


RT 


x 2 (l - x 2 ) 


- 2 RTA 


(7.5) 
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The system is unstable when 


A > 


1 

2x 2 (1 - X 2 ) 


(7.6) 


The right-hand side of Eq. (7.6) is plotted in Fig. 7.3. In Fig. 7.3, the x-axis represents the mole 
fraction of component 2, and the y - axis represents the value of the parameter A in the Margules model. 
The shaded region denotes the area where the inequality given in Eq. (7.6) is true, and therefore, the 
system will undergo liquid-liquid phase separation. Outside this region, the mixture is stable. 
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Figure 7.3: Stability diagram for the two-suffix Margules equation. The shaded region represents 
states where the Margules model predicts that the system is unstable. 


7.2 Liquid-liquid equilibrium 


When a liquid mixture becomes unstable, it separates into two coexisting liquid phases. We label one 
of the phases with a ' and the other phase with a // . The criterion for phase equilibrium is 


M a Ma 

Ma + RTln x'ala = Ma + RT l rl x 'all 
RT In x ^' 0 = RT In®"- 7" 

Y y = y'V' 

/a /a 


(7.7) 


Given a model for the activity coefficients, Eq. (7.7) yields a set of constraints that can be solved for 
the compositions of each liquid phase. 

Conversely, given the composition of the coexisting liquid phases, Eq. (7.7) can also be used to fit an 
activity coefficient model. For example, using the two-suffix Margules model for a binary mixture, 
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we can determine the parameter A 


ln -^77 = in — 7 

ii % 


Ax'2 — Ax'2 2 = In 


1 

x'{ 

/ 

1 


X 


A = 


x ' 2 — x\ 


t" 

In ^ 

2 x l 


Jt 2 


( 7 . 8 ) 
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8 Solid-liquid equilibrium 

8.1 Introduction 

At atmospheric pressure, water will freeze at 0°C. If a solute is added to water, then the mixture will 
freeze at a lower temperature. This is known as freezing point depression and occurs with any general 
mixture. Using a model for the activity coefficients, we can quantitatively predict the magnitude of 
this effect. 


8.2 Phase behavior 

To demonstrate the general freezing behavior of a binary mixture, we present the solid-liquid phase 
diagram for systems of ethanol and water at 1 atm in Fig. 8.1. The solid line with the filled symbols 
is the freezing curve of water in the mixture. Above the curve the solution is completely liquid; 
below the curve, it is a liquid mixture coexists with solid water (i.e., ice). At this pressure, pure water 
freezes at 273.15 K. As ethanol is added to the solution, the temperature at which ice begins to form 
in gradually decreases. 

The dotted line with the open symbols on the right side of the diagram is the freezing curve of 
ethanol in the mixture. As one passes downward through this curve, solid ethanol precipitates from 
the solution. Note that the temperature at which ethanol begins to freeze decreases as the amount 
of water in the solution increases. Therefore, we see that the addition of impurities decreases the 
freezing temperature of a substance. This phenomena is known as freezing point depression. 

The point where the two freezing curves intersect with each other is called the eutectic point. This is 
the lowest temperature at which a binary mixture can remain in liquid without either of its components 
precipitating. 

Finally, we note that below the dotted line in Fig. 8. 1 the system exists as two coexisting solid phases: 
one consisting of pure ice, and the other composed of pure solid ethanol. The relative amounts of 
these two phases is given by the lever rule. 


8.3 Conditions for equilibrium 

In this section, we derive the equation that governs the shape of the freezing curve. Let’s consider the 
freezing of a species in a general multicomponent mixture. We make the assumption that the solid 
phase consists of pure component a. At equilibrium, the chemical potential of the solid phase is the 
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Figure 8.1: Solid-liquid phase diagram for mixtures of ethanol and water at 1 atm. Data taken from 
SU Pickering, J. Chem. Soc. 63, 998 (1893). 


same as that in the liquid phase: 


u (9 = M 


+ RT\nx ala = M°’ (s) (T,p) 

RT\nx ala = M°’ (s) (T,p) - lx°J l \T,p) 
AG a (T,p) 


lnx a 7 a = 


RT 


( 8 . 1 ) 


where A G a is the difference in the chemical potentials of the solid and liquid phases (i.e., A G a (T, p) = 
li « ( s ) (T, p) — fia^ (T, p)). If we knew the composition of the mixture (i.e., x a ), 7 a , and the function 
A G a , then we could use Eq. (8.1) to predict the freezing temperature T of the mixture. 


One method to approximate the function A G a is by using a Taylor series expansion around the 
temperature T m , the melting temperature of pure a at pressure p. This yields 

A G a {T,p) « A G a (T m ,p) + dAGa ^ m,P \ T - T m ) + • • • (8.2) 


At the melting temperature T m , the molar Gibbs free energy of the pure solid phase is equal to the 
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molar Gibbs free energy of the pure liquid phase. That is 

AG a (T m ,p) = 0 


(8.3) 


Therefore, the first term of the Taylor series expansion in Eq. (8.2) vanishes. 


From the fundamental equation of thermodynamics, we arrive at 


dAG a 

dT 


-A S a 


(8.4) 


where A S a is the molar entropy change of freezing. From the relation G = H — TS , we can relate 
the entropy change of freezing to the enthalpy of freezing A H a , which is the heat released by the 
system when it transforms from a liquid to a solid: 


AG a (T m ,p) = A H a (T m ,p) - T m AS a (T m ,p) 
0 = A H a (T m ,p) - T m AS a (T m ,p) 

AS a (T m ,p) = AH a (T m ,p ). 

m 

Combining this with Eq. (8.2) yields 


A G a (T,p) ~ 


A H a (T m ,p) 


(T - T m ) + 


(8.5) 


( 8 . 6 ) 
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Finally, we find: 


In XoTIcl = 
In XoTIa = 


AH a (T m ,p) 


RTT m 


(T - T m ) 


A H a (T m ,p) 


RT„ 



(8.7) 


This equation allows the prediction of the freezing curve of a mixture. It requires from knowledge of 
the freezing temperature and enthalpy of melting of the pure component a, as well as a model of the 
activity coefficients of the liquid mixture. In the absence of information for the activity coefficients, 
the ideal solution model can be used (i.e., = 1); Eq. (8.7) then reduces to the van’t Hoff equation. 
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9 Gas solubility and Henry’s law 


In order to describe vapor-liquid equilibria, we have been using Raoult’s law (see Eq. (6.1)) or the 
modified Raoult’s law (see Eq. (6.3)). One of the required inputs to the modified Raoult’s law is the 
vapor pressure p™ p (T) of each component a in the pure state; this is the pressure below which pure 
liquid a will begin to boil (or vaporize) at a given temperature T. However, the vapor pressure of a 
substance is only defined at temperatures below its critical temperature. Above the critical temper- 
ature, the vapor pressure is undefined (recall that the vapor pressure curve terminates at the critical 
point). Therefore, Raoult’s law and the modified Raoult’s law are only applicable if the critical tem- 
peratures of all the components in the system are above the system temperature. 

This presents difficulties when we want to examine the solubility of gases, such as oxygen which has a 
critical temperature of T c = 154.59 K and nitrogen which has a critical temperature of T c = 126.21 K, 
in liquids. The critical temperature of these gases are well below ambient temperatures. For these 
systems, another approach is required to determine the conditions for vapor-liquid equilibrium. 


9.1 Henry’s law 

Let’s consider a binary mixture of carbon dioxide and hexane at a temperature of 393.15 K. The 
critical temperature of hexane is 507.5 K, which is higher than the system temperature, so pure liquid 
hexane can exist as a liquid. The critical temperature of carbon dioxide is 304.2 K, which is lower than 
the system temperature. Consequently, pure carbon dioxide does not possess a vapor-liquid transition 
at this temperature, and its vapor pressure is undefined. Raoult’s law cannot be applied to this system. 

In Fig. 9.1, we plot the variation of the partial pressure of carbon dioxide with the mole fraction of 
the gas dissolved in a liquid mixture of carbon dioxide and hexane. When xco 2 = 0, the system 
consists solely of hexane, and the partial pressure of carbon dioxide is zero. The vapor space consists 
entirely of hexane at its vapor pressure. As the concentration of carbon dioxide dissolved in the liquid 
phase increases, the partial pressure of carbon dioxide in the vapor phase also increases. However, 
for xqo 2 > 0.7, the liquid phase no longer exists; above this mole fraction, the system is a one phase 
vapor (or gas) mixture. 

For this system, we cannot apply Raoult’s law, as the vapor pressure of carbon dioxide is undefined 
because the system temperature is above its critical temperature. To overcome this difficulty, for each 
component that has a critical temperature below the system temperature, we replace Raoult’s law with 
Henry’s law, first proposed by J. W. Henry in 1800 to empirically describe the solubility of gases: 

p a = K a (T)x a , (9.1) 

where K a ( T ) is known as the Henry’s law constant. 

The Henry’s law constant is equal to the slope of the partial pressure of component a at infinite 
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Figure 9.1: Variation of the partial pressure of component carbon dioxide in the vapor phase with 
its mole fraction in the coexisting liquid phase mixture of carbon dioxide and hexane at 393.15 K. 
The solid line and filled symbols denote the actual variation of the partial pressure. The dotted line 
is Henry’s law. Data taken from YH Li, KH Dillard, and RL Robinson, J. Chem. Eng. Data. 26, 53 
(1981). 


dilution. Therefore, Henry’s law is expected to be accurate where the concentration of a in the liquid 
mixture is low. K a (T ) plays analogous role as the vapor pressure in Raoult’s law (see Eq. (6.1)). 


9.2 Activity coefficients 


As with Raoult’s law, Henry’s law is only an approximation. However, we can generalize its range of 
validity by defining an activity coefficient 7 * as 


* 

7a 

* 

7a 


actual p a 
ideal p a 

Pa 


(9.2) 

(9.3) 


Note that Henry’s law becomes exact in the limit x a — ► 0 and, therefore, we expect that 7 * — ► 1 as 
x a — ► 0. This is a different normalization from 
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Figure 9.2: Variation of the partial pressure of component a in the vapor phase with its mole fraction 
in the coexisting liquid phase. The solid line is the actual variation of the partial pressure. The dashed 
line is the prediction of Raoult’s law, and the dotted line is the prediction of Henry’s law. 


In Fig. 9.2, we plot a schematic drawing of the variation of the partial pressure of a component with 
its mole fraction in the liquid phase. The solid line represent the experimental variation of the partial 
pressure, the dashed line is the prediction of Raoult’s law, and the dotted line is the prediction of 
Henry’s law. 

Even if a component has a vapor pressure at the system temperature (i.e., its critical temperature is 
greater than the system temperature), we still may prefer to use Henry’s law instead of Raoult’s law. 
This is usually done for components that are very dilute, where Henry’s law becomes more and more 
accurate. 

If we compare Henry’s law and the modified Raoult’s law, then we find that Henry’s law can be 
considered as the modified Raoult’s law with 7 ^ = K a /p™ p . 
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10 Equations of state 


An equation of state (EOS) provides a relation between the pressure, temperature, volume, and com- 
position of a system. Typically, this is given as an explicit function of temperature, molar volume, 
and composition p = p(T, V, X 2 , ^ 3 , • • • ). One example is the ideal gas law, where the pressure is 
given by 

TDJi 

p=— ( 10 . 1 ) 

where R is the universal gas constant. 

Real materials have much more complicated equations of state. In this chapter, we present other 
equations of state that are used to model real systems and demonstrate how to use them to determine 
various thermodynamics properties from them. 


10.1 The principle of corresponding states 


The properties of different fluids are in general quite different from each other. In Fig. 10.1, we plot 
the liquid-vapor coexistence curves for a variety of fluids. However, Guggenheim observed that if 
the properties of a certain class fluids are plotted in terms of reduced variables, variables divided by 
their critical values, the data can be collapsed onto a single curve, independent of the identity of the 
fluid. This can be seen on the right plots in Fig. 10.1. This collapse of the properties is known as the 
principle of corresponding states. 


This observation implies that the thermodynamic properties of all materials are given by a single, 
“universal” function, independent of the identity of the material. For example, defining the compress- 
ibility factor Z as 


Z = 


pV 

Sr’ 


the principle of corresponding states implies that 


( 10 . 2 ) 


Z = Z(T r ,p r ) (10.3) 

where T r = T/T c is the reduced temperature, p r = p/p c is the reduced pressure, T c is the critical 
temperature of the material, and p c is its critical pressure. 


Guggenheim observed that this principle actually holds very well for many simple fluids, fluids com- 
posed of molecules that are roughly spherical. However, in general, this principle does not hold for 
an arbitrary fluid, in particular, for polar fluids, such as water, and for fluids composed of highly 
non-spherical molecules, such has n-hexane. 


Pitzer noticed that as the molecules in a fluid begin to elongate and deviate from a spherical shape, 
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Figure 10.1: Vapor-liquid coexistence curves for neon (open squares), argon (filled circles), 
krypton (filled squares), xenon (open down-triangles), methane (open circles), nitrogen (open 
diamonds), and oxygen (open up-triangles). Data taken from the NIST Chemistry Webbook 
http : //webbook . nist . gov. 


the slope of the vapor pressure curve increased in magnitude. In order to quantify this effect, Pitzer 
defined the accentric factor u as 


u = -1 -log 10 


p va P ( T = 072;) 

Pc 


(10.4) 


For spherical molecules, u is approximately equal to zero. 


This generalized principle of corresponding states implies that the compressibility factor of a fluid 
can be expressed as 

Z — Z(T r ,p r ,ou) (10.5) 

Various approximations for this function have been developed in the literature and are presented 
in standard reference texts. Therefore, given this function, we can determine the thermodynamic 
behavior of a compound from knowledge of its critical temperature, critical pressure. 
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10.2 The van der Waals equation and cubic equations of state 


The basic assumption in deriving the ideal gas equation of state is that the interactions between the 
molecules that compose the gas are negligible. At very low densities, where the average distance 
between molecules is extremely large compared to the range of the intermolecular interactions (or 
size of a molecule) and collisions between molecules occur only rarely, this assumption is fairly good. 
However, for moderately dense gases, this assumption begins to breakdown, as collisions between 
molecules become more and more frequent. In the case of liquids, this assumption is completely 
erroneous, as intermolecular interactions dominate the properties of the liquid. 


To account for the finite size of molecules in the system. If we consider molecules to be of volume 6, 
then the actual volume available to the system is approximately V — b. Therefore, we replace V in 
the ideal gas equation of state with V — b. In addition, at large distances, the molecules in a gas exert 
an attractive interaction between each other due to dispersion forces. The van der Waals equation 
of state, which accounts for both the short-ranged repulsions and long-ranged attractions between 
molecules, is given by 


RT a 
V -b ~ V 2 


( 10 . 6 ) 


A schematic plot of the variation of the pressure with volume, as predicted by the van der Waals 
equation of state, at various temperatures is given in Fig. 10.2. At temperatures above the critical 
temperature, the pressure-volume variation is monotonic and qualitatively similar to that of an ideal 
gas (see dotted-line). At temperatures below the critical temperature, the pressure-volume curve be- 
gins to oscillate, exhibiting a “van der Waals” loop (see dashed-line). This behavior is unphysical, but 
represents the vapor-liquid transition, and should be replaced by the solid line. The precise location 
of the solid line is given by the Maxwell construction. 


When the temperature is equal to the critical temperature, the pressure-volume curve exhibits a point 
with the inflection (concavity) of the curve changes sign. This occurs at the critical point, where we 
have: 



d 2 p 

W 2 


= 0 


(10.7) 

( 10 . 8 ) 


Applying Eqs. (10.7) and (10.8) to the van der Waals equation of state, we find 

dp = RT C 2 a 

dV ~ (V c - b) 2 y c 3 
d 2 p _ 2 RT 6a _ 

W* ~ (v c - by ~ vy ~ 0 


(10.9) 

( 10 . 10 ) 
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Figure 10.2: Variation of pressure with volume for the van der Waals equation of state for: (i) T >T C 
(dotted line), (ii) T — T c (solid line), and (iii) T <T C (dashed line). 
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Multiplying the Eq. (10.9) by 2/ (V c — 6) and adding it Eq. (10.10) yields 
2 RT 4 a 2 RT 6 a _ _ 

~Wc - w + mvc z r + Wc 

4 a 6 a n 

mVc-b)~V^ = 

2V C = 3(F C - b) 
v c = 3b 


Substituting V c = 3 b into Eq. (10.9) yields 


RT C ( 2 a 

(Vc - W + V? 

nr. 


0 

2 a(V c - b) 2 

K 3 

2a(36 — b ) 2 

8 a 
27 b 


( 10 . 11 ) 


We can compute the critical pressure by substituting T c and V c into the equation of state, 

RT a 
Pc = V^b ~ V 2 
8a/(276) a 

= 3b -b ~ (3 W 
a 

~ 2762 


( 10 . 12 ) 


We have shown that, given the parameters a and b of the van der Waals equation of state, one can 
determine the critical properties of a fluid. However, we can perform the inverse task. That is, given 
the critical temperature and pressure for a fluid, one can determine the a and b parameters of the van 
der Waals equation of state. Dividing Eq. (10.1 1) by Eq. (10.12) 


RT C _ 8a/ (276) 
Pc a/ {27b 2 ) 

7 RT C 
b = — — 

8p c 

Inserting the expression for 6 into Eq. (10.1 1), we find 


27RT c /{8p c ) 
27 (RTcf 
64 p c 


(10.13) 


(10.14) 

(10.15) 


74 


Download free eBooks at bookboon.com 


Chemical Thermodynamics 


Equations of state 


Therefore, given the critical temperature and pressure of a fluid, we can use Eqs. (10.13) and (10.15) 
to determine the parameters for the van der Waals equation. This then allows us to make predictions 
about the thermodynamic behavior of the fluid at any other state point. 

Many other equations of state exist beyond the van der Waals EOS. Some examples include: 


• Redlich-Kwong equation of state 


RT a 

P ~ V - b ~ T^W(y + b ) 

• Soave-Redlich-Kwong equation of state 

RT a(T) 

P ~ V-b ~ F(F + b) 


where 


a(T ) = a(T c )a(T) 
a{T) = [1 + ft(l - r^ 2 )] 2 

ft = 0.480 + 1.574a; - 0.176a; 2 


(10.16) 


(10.17) 


(10.18) 


• Peng-Robinson equation of state 


where 


RT a(T) 

V-b ~ F 2 + 2 W - 6 2 


a(T) = a(T c )a(T) 

a(T) = [1 + k( 1 - T r 1/2 )] 2 

ft = 0.37464 + 1.54226a; - 0.26992a; 2 


(10.19) 


The relationship between the parameters of these models and the critical temperature and pressure of 
the fluid are summarized in Table 1. 


Table 1 : Parameters for various equations of state in terms of the critical temperature and pressure. 


a(T c ) 

van der Waals i? (RTc) 

64 p c 

Redlich-Kwong 0.42748 

Redlich-Kwong-Soave 0.42748 
Peng-Robinson 0.45724 


b 

RTc 

8p c 

0.08664^ 

Pc 

0.08664 4^ 

Pc 

0.07780 ^ 

Pc 
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10.3 Equations of state for mixtures 

Thus far, we have only been considering single component systems. One method of extending an 
equation of state to mixtures is by making the parameters a and b a function of composition. This 
is done through the application of mixing rules. For example, we can define the parameter b of the 
mixture as 


b = Y , y « h « 00.20) 

OL 

where b a is the b parameter for pure a , and y a is the mole fraction of a in the system. 

Similarly, we can define the parameter a of the mixture as 

a = ^ y « y “ ,cw (10.21) 

a, a' 


where a aa t = ( 1 — k aa r ) (a a a a / ) l / 2 , a a is the a parameter of pure a , and k aa * is the binary interaction 
parameter between species a and a'. The binary interaction parameter k aa r is a purely empirical 
parameter that is used to obtain better agreement between the predictions of the equation of state and 
experimental data. Typical values of the parameter are given in Table 2. 

Table 2: Typical values of the binary interaction parameters for various mixtures. 


mixture 

to' 

hydrocarbon/hydrocarbon 

~ 0 

CCVhydrocarbon 

~ 0.15 

water/hydrocarbon 

~ 0.5 
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11 Thermodynamics from equations of state 


Equations of state relate the pressure, temperature, volume, and composition of a system to each 
other. In this Chapter, we show how to determine other thermodynamic properties of the system 
from an equation of state. In a typical equation of state, the pressure is given as an explicit function 
of temperature, volume, and composition. Therefore, the “natural” variables are the temperature, 
volume, and composition of the system. That is, once given the volume, temperature, and composition 
of the system, the pressure is readily calculated from the equation of state. 

Once one of the free energies of a system is known as a function of its natural variables, then all 
the other thermodynamic properties of the system can be derived. For these equations of states, the 
Helmholtz free energy is the relevant quantity. In the following, we demonstrate how to determine 
the Helmholtz free energy from an equation and then proceed to show how to derive other properties 
from it. 


11.1 The residual Helmholtz free energy 

The free energy that has temperature, volume, and mole numbers as its natural variables is the 
Helmholtz free energy. Before we stated that once the Gibb’s free energy of a system is known as a 
function of temperature, pressure, and mole numbers G(T,p, TVi, N 2 , . . . ), all the thermodynamics 
of the system are known. This is equivalent to the statement that once the Helmholtz free energy is 

known as a function of temperature, volume, and mole numbers of the system A(T, V, TVi, ) ? 

all the thermodynamics of the system are known. The fundamental equation of thermodynamics can 
be written in terms of the Helmholtz free energy as 

dA = —SdT — pdV + ^ Ha dN a (11.1) 

OL 


Given and equation of state, we can determine an explicit expression for the Helmholtz free energy 
from the fundamental equation of thermodynamics. At constant temperature and mole numbers, we 
have 


dA — —pdY 


( 11 . 2 ) 


By integrating both sides of this equation from a total volume Yi to a total volume Y 2 

r ¥. 2 
'¥1 


r ¥2 

A(T,Y2,N u N 29 ...)-A(T,Yi,N 1 ,N 2 ,...) = - / dY'p 

Jv 1 


(11.3) 
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For an ideal gas, we have p = NRT/V , so 


A l9 (T, V 2 , N 1} N 2 ,...)~ A i9 (T, Y u N u N2, ...) = 


rY, 

JVx 


5/2 ,NRT 


dV' 


¥' 


(H.4) 


NRT\ 

“T 77 " ) 


Subtracting Eq. (1 1.4) from Eq. (1 1.3), we find 

[A(T, V 2 , N\, N 2 , . . . ) — A(T, Yi,N 1 ,N 2 , . . . )] 

~[A i9 (T, y 2 , N U N 2 , ...)- A i9 (T, Vi,N 1 ,N 2 , ...)] = - jj 2 dV ^ ■ 

[A(T, Y 2 , N\ , N 2 , . . . ) A i9 (T , ¥ 2 , N u N 2 , . . . )] 

-IA(T, ¥i,Ni,N 2 , ...)- A i9 (T, Yi,Ni,N 2 , ...)] = -NRT ^ ^{Z - 1) (1 1.5) 

A res (T, ¥ 2 , N u N 2 , A res (T, Yi,N u N 2 , ...) = -NRT R %-[Z - 1) 

JYi ¥ 

where we have defined a residual property X res as the difference of the value of the property of the 
system and that of an ideal gas at the same temperature, volume, and species mole numbers 


X res {T,V,N 1 ,N 2 ,...) = X(T,V,N 1 ,N 2 ,...)-X i9 (T,V,N 1 ,N 2 ,...) (11.6) 


As the density of the system approaches zero, or the molar volume of the system approaches infinity, 
the properties of a system approach those of an ideal gas. In particular, 

A(T,¥,N 1 ,N 2 ,...)^A i9 (T,¥,N 1 ,N 2 ,...) as ¥ - oo. (11.7) 

If we let ¥2 = Y and \ 1 — x in Eq. (1 1.6), then we find 

f V dv' 

A res (T,¥N 1 ,N 2 ,...) = ~NRT -fr( z - 1) 

Toe ¥ 

dV f 

= NRT Jv Y r{z ~ 1) (1L8) 

Given an expression for the compressibility factor Z as a function of temperature and volume for a 
fluid (an equation of state), we can determine the residual Helmholtz free energy of the system. 


78 


Download free eBooks at bookboon.com 


Chemical Thermodynamics 


Thermodynamics from equations of state 


For the van der Waals equation of state, the compressibility factor is given by 

pV 


z = 


RT 

V ( RT 


Z- 1 = 


RT VV 

1 

cr 1 

1 

'SI 

V 

a 

1 

RTV 

b 

a 

1 

RTV 


(11.9) 


Substituting this expression into Eq. (1 1.8), we can obtain the residual Helmholtz free energy for the 
van der Waals equation of state: 


NRT 


A res (T, V, N u N 2 , . . . ) _ f°° dV 

~ Jv 'V 

= / 

Jv 


dV 


V-b RTV / 

1 1 a 


\ V-b 
In — — 1- 


V - b V RTV 2 

OO 


(11.10) 


a 


= In 


V 

V 

V-b 


RTV 

a 


v 


RTV 


The residual chemical potential of species a can be obtained from the residual Helmholtz free energy 
by differentiating by the number of moles of a (see Eq. (11.1)): 




res 

a 


dA res 

~dN ~ 


( 11 . 11 ) 


For the van der Waals equation of state, this becomes 


res 

ra 


d 


dN a 
= RT In 


NRT 

V 


In 


V 


V~Nb 


Na ' 
RTV 


V-Nb 


+ 


NRT \ dNb 


V-Nb dN, 


1 dN 2 a 
V dN^ 


( 11 . 12 ) 


To obtain the derivative of the a and b parameters with respect to N a , we need to use the mixing rules 
given in Eqs. (10.20) and (10.21). This yields 


dNb 

dN~ a 


d 

d 

dN a 


N ^ ^ Ua'ba' 
a' 

^ ^ Nql' b a r 
a ' 


(11.13) 
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and 


dN 2 a 


d 


dN a 

d 

dN n 


JV 2 53 


Ua'lla'^a'a" 


J2 N « iN < 


ol" 


2j>. 


7 CLcvry ' 


2 


ya'dc 


( 11 . 14 ) 


Using these expressions, we find that the residual chemical potnetial for the van der Waals equation 
of state is given by 


y 


res 

a 


RT In 


V 


V-b 


RTb a 

-\ 

V-b 



( 11 . 15 ) 
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11.2 Fugacity 


An ideal gas is an ideal solution. Thus, the chemical potential of species a in an ideal gas mixture at 
temperature T, pressure p, and species mole fractions y\, . . is given by 


Ma(T>.P> 2 /l, 2 / 2 , • • • ) = Ma lg ( r >p) + RT^Va 


where pa(T,p) is the chemical potential of pure a. 


(11.16) 


The difference in the chemical potential of an ideal gas at pressure p and mole fractions yi,y2, ••• and 
an ideal gas at pressure p' and mole fractions y[,y r 2 ,... is given by 


t%(T,P, 2/1, 2/2, . . . ) = /£ KT,p\ y'l, 2/2, • • • ) + RT In ^ 

= ^ S T,p',2/i,2/2, . . . ) + RTln^ (11.17) 

Pa 

Note that we have made use of the relation In analogy with ideal gases, we define the fugacity of a 
real gas as 


Va(T,p, 2/1, 2/2, • • • ) = Va(T,p', y[, 2/2, • • • ) + RT In 


/q(T,p, 2 / 1 , 2/2, •••) 
fa(T,p', y[,y' 2 , ■ ■ •) 


(11.18) 


where f a is a function of temperature, pressure, and the system composition. 


The fugacity coefficient <f a is defined as 


0a 


fa 

Pa 


fa 

VaP 


(11.19) 


This gives us the difference between the chemical potential of the system and the chemical potential 
of an ideal gas at the same temperature, pressure, and composition. That is: 

y a (T,p,y 2 ,V 3 , ■ ■ ■) - l^a(T,p,y 2 ,y 3 , ■ ■ ■) = RT\n.(j) a (11.20) 

The residual chemical potential is equal to the difference between the chemical potential of the system 
and the chemical potential of an ideal gas at the same temperature, molar volume, and composition. 
The difference given in Eq. (1 1.20) is between the system and an ideal gas at the same temperature, 
pressure , and composition. This is not the same as the residual chemical potential. However, the 
two quantities are related. If we consider an ideal gas with the same molar volume as the system, its 
pressure will be equal to p' = RT /V , which will, in general, not be equal to the actual pressure p of 
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the system. 


fJ>a(T,p,y 2 ,y 3 ,...) = H l £(T,p' = RT /V, y 2 , y 3 , ■ ■ ■ ) + RTIn 

= ^(T,p' = RT/V,y 2 ,y 3 ,.. 

= /#(' T,p ' = RT/V , y 2 , y 3 ,...) + RTlnZ 


pV 

+ RT ln^— 

RT 


( 11 . 21 ) 

(11.22) 


where Z is the compressibility factor of the system. Substituting this expression in to Eq. (1 1.20), we 
find: 


RT In 4> a = Va(T,p,y 2 ,y 3 ,...) ~ Va(T,p,y 2 ,y3,...) 

= Va{T,p,y 2 ,y 3 ,...) - y^(T,p',y 2 ,y 3 ,...) - RT In Z 
= ^ s (T,p,y 2 ,y 3 , RT In Z (11.23) 

Therefore, given an equation of state, we can use Eq. (1 1.23) to obtain an expression for the fugacity 
coefficient of a system. 


11.3 Vapor-liquid equilibrium with a non-ideal vapor phase 

Raoult’s law assumes that the liquid phase is an ideal solution and the vapor phase is an ideal gas. 
With an equation of state, such as the van der Waals equation of state, the vapor phase can now 
described more realistically, and we do not need to make the ideal gas assumption. In this situation, 
the chemical potential of the vapor phase 


»,...) = PL% C7>™ p ) + Mr, M2, . . . ) - K(r,K* p ) 

J ax- 1 - ? Pol ) 

where /° is the fugacity of pure component a. 

Now if we equate the chemical potentials of the liquid and vapor phases, we find 


(11.24) 


y 9 a (T,p,y 2 ,...) = fi l a (T,p, x 2 , . . . ) 
l-^ 9 (TlC P ) + RT\n f f o = ii°J(T,p) + RTln la x a 

J a 


f C 

(C s (r,p“ p ) + RT In f o = PJ(T, P ™V) + / 

J a J p 


d P O 

p dp 


+ RT In 


RT In = 


RT In 


y a p(f>a(T,p,y 2 ,...] 

P rr a (T, P v a ap ) 


f Pa i 

/ dp'V°’ l (T,p) + RT ln la x a 
Jp 

fp V T P 

/ dp'V°\T,p)+RT\n la x a 
Jp 
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Neglecting the first term on the right side of the above equation is equivalent to assuming that the 
Poynting factor is equal to one. Making this assumption leads to 


K ap ^(T,K ap ) 

<i>cx{T,p, J/2) • • • ) 

VaP 0°(T,K ap ) 




laX a P 


vap 

a 


(11.25) 


If the vapor phase is ideal, then the fugacity coefficients are equal to one, and we see that Eq. (1 1.25) 
becomes the modified Raoult’s law that we have derived previously. 


We do not need to apply an activity coefficient for the liquid phase. In fact, we can, in principle, apply 
the equation of state to determine the thermodynamics of both phases. In this case, we have 


Ma = l& 


Va(T,P, V 


+ RT In 4> 9 a = fj^(T,p,xi,.. 
Vo.pK = x a p<f) l a 


+ RT In < 


(11.26) 


where we have used Eq. (11.17). Finally, we note that the ratio of the vapor mole fraction to the liquid 
mole fraction is related to the ratio of the fugacity coefficients: 


Va <t> l a 



(11.27) 
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12 Chemical reaction equilibria 


In this section, we discuss the thermodynamics of systems undergoing chemical reactions. First we 
derive the conditions for equilibrium. Then, we generalize the Gibbs phase rule for reacting systems. 
Afterwards, we examine the calculation of the and introduce the equilibrium constant. Finally, we 
examine the influence of pressure and temperature on chemical reaction equilibria. 


12.1 Conditions for equilibrium 

Let’s consider a system that initially contains moles of nitrogen, moles of hydrogen, and 
^n°h 3 m °l es °f ammonia. This system can undergo the following reversible chemical reaction: 


N 2 + 3H 2 ^2NH 3 (12.1) 

where nitrogen and hydrogen react to form ammonia, or ammonia can decompose to form nitrogen 
and hydrogen. 

If the reaction proceeds forward by £ moles (i.e., £ moles of nitrogen combine with 3£ moles of 
hydrogen combine to form 2£ moles of ammonia), then, from a mole balance, the number of moles 
of each species in the system is given by 

N. n 2 = < 2 } -£ 

Nh 2 = - 3 £ ( 12 . 2 ) 

N NK3 = N$l 3 + 2 £ 

The quantity £ is known as the extent of reaction. If the system does not react, then the extent 
of reaction equals zero (i.e., £ = 0), and the mole numbers are equal to their initial value (i.e., 
N a = N^)- If the reaction proceeds backwards (i.e., ammonia decomposes into nitrogen and 
hydrogen), then £ has a negative value. 

Note that the coefficient in front of the extent of reaction in the second term on the right side of 
Eqs. (12.2) is known as the stoichiometric coefficient for the component in the reaction. For a general 
reaction, we have 


N a = (12.3) 

where u a is the stoichiometric coefficient for component a in the reaction. The stoichiometric coef- 
ficient is positive for products of the reaction, and it is negative for the reactants of the reaction. The 
absolute value of the stoichiometric coefficient of a component is equal to the numerical value of its 
coefficient in the reaction. For example, in the reaction given in Eq. (12.1), we have 

= -1, ^h 2 = z/nh 3 = 2 (12.4) 
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From Eq. (12.3), we directly relate the changes of the number of moles of species a with the change 
in the extent of reaction: 

dN a = v^di (12.5) 

For an isolated system where the total energy and volume are constant (i.e., dU = 0 and Y = 0), the 
total entropy is maximized at equilibrium; therefore, we must have dS = 0. The change in entropy 
of a system with respect to changes in the energy, volume, and number of moles of each species is 
given by (see Eq. (3.1)) 

dS=^dU + ^dY-J2Y dNa 

a 

a 

where we have used Eq. (12.5) to relate changes in the mole numbers to the extent of reaction. 

In order to ensure that dS = 0, the coefficient of the d £ term must equal zero. Consequently, the 
criterion for chemical reaction equilibrium is 

y% a // a = 0 (12.7) 

a 

for each chemical reaction that takes place in the system. 
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12.2 The phase rule for chemically reacting systems 

Because the presence of chemical reactions introduces additional constraints on a system, we need 
to reconsider the Gibbs phase rule. In the absence of chemical reactions, the number of degrees of 
freedom / in a system is given by 


/ = 2 + U — 7T 

where uo is the number of components, and n is the number of phases. 

For a system with chemical reactions, there is one addition equation (see Eq. (12.7)) for every in- 
dependent chemical reaction in the system. The presence of chemical reactions, therefore, imposes 
constraints on a system, which decrease its degrees of freedom. The phase rule now becomes: 

f — 2 -\- (jo — tv — r (12.8) 

where r is the number of independent reactions. 


12.3 Gas phase reactions 


Now we discuss how to calculate the equilibrium composition of a reacting non-ideal gas mixture. In 
order to use Eq. (12.7), we need an expression for the dependence of the chemical potential of each 
species in the system with its composition. In general, this can be expressed as (see Sec. 1 1): 

Va(T,p,y 2 ,y 3 ,.. . ) = /#- °(T,p 0 ) + RT \ n y J^P (12.9) 

Po 

where fia°(T^p) is the chemical potential of pure a in an ideal gas state at temperature T and 
pressure p, and po is a reference pressure (typically taken to be either 1 atm or 1 bar). Inserting this 
expression for the chemical potentials into the condition for chemical reaction equilibrium Eq. (12.7), 
we find 


^ ^ L'afJ'a — 0 
a 

Va[p$'°(T,p) + RT In y a cf> a p} = 0 

a 

RTlnl[(y a ct> a pr «* = - ^ a /# 0 (T,p) 

a a 

f \{y a <t>aP) Va = K eq (T,p 0 ) 

a 


( 12 . 10 ) 
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where K eq is the equilibrium constant, defined as: 


Keq(T,p 0 ) = exp 


AG rxn (r,p 0 ) 

RT 


and AG rxn (T,po) = Y^ a u atfa ,0 (T,pQ) is the Gibbs free energy of reaction. 


( 12 . 11 ) 


The equilibrium constant and the Gibbs free energy of reaction are independent of the composition 
and pressure of the system but are dependent on the temperature of the system and the choice of the 
reference pressure. The greater the value of the equilibrium constant (which corresponds to more 
negative values of the Gibbs free energy of reaction), the further the reaction proceeds to completion. 
Given a value for the equilibrium constant and the pressure of the system, Eq. (12.10) can be solved 
to determine the equilibrium composition of the system. In the next section, we discuss how to 
determine the value of the equilibrium constant. 


12.4 The standard Gibbs free energy of formation 

In order to determine the equilibrium constant, we need to be able to compute the Gibbs free energy 
of reaction AG rxn * This is dependent on the quantities fia ’°(T,po) — the chemical potential of pure 
species a in an ideal gas state at temperature T and pressure po- These can be directly related to the 
standard Gibbs free energy of formation. 

The standard Gibbs free energy of formation of a compound is defined as the Gibbs free energy 
change for the reaction of forming one mole of the compound from its constituent elements in their 
standard state at a temperature of 25° C and a pressure of 1 bar (or 1 atm). The standard state of an 
element is its stable form at 25°C and 1 bar. For example, the standard state for oxygen is the form 
O 2 in the (ideal) gas state, the standard state for bromine is in the liquid state, and the standard state 
for carbon is in the form of graphite. The standard Gibbs free energy of formation of elements in their 
standard state at 25°C and 1 bar is defined to be zero. Tabulations of values for the standard Gibbs 
free energy of formation for a broad range of compounds can be found in standard references, such 
as the CRC Handbook of Chemistry and Physics or Perry’s Handbook of Chemical Engineering, or 
on online resource, such as the NIST Chemistry webbook (http : //webbook . nist . gov). The 
standard enthalpy and Gibbs free energy of formation are listed in Table 3. 

Table 3: Standard formation enthalpy and Gibbs free energies of various compounds at T = 25° C 
and p— 1 bar. 



Hf 

Gf 


(kJmol 1 ) 

(kJmol 1 ) 

H 2(0) 

0.00 

0.00 

N2G7) 

0.00 

0.00 

NH 3 ( 0 ) 

-46.11 

-16.45 
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The Gibbs free energy of the reaction in Eq. (12.1) can be obtained from Table 3 

AG rxn (25°C,lbar) = 2(-16.45) - 1(0.00) - 3(0.00) = -32.9 kJ mol" 1 , 

From this, we can compute the corresponding value of the equilibrium constant. 

Tables of the standard Gibbs free energy of formation allow the calculation of the equilibrium constant 
at standard conditions (i.e., T = 25°C and p = lbar). However, in most situations, we require 
the equilibrium constant at other conditions. In the next section, we discuss how to estimate the 
equilibrium constant at other temperatures. 


12.5 The influence of temperature 


Often, we want to know how the equilibrium constant depends on temperature to determine the opti- 
mal conditions to perform a reaction. To determine this, we need to know how AG rxn / (RT) changes 
with temperature. The variation of the Gibbs free energy with temperature is given by 


which leads to 


(G\_dG GdT 

\t) ~ ~t - ~T r 

_ —SdT + Vdp (H - TS)dT 

J 1 ^2 

H V 
= Hd (2) + ^ dp 


d(G/T) 

9(1 IT) 


( 12 . 12 ) 


(12.13) 


Therefore, we find that the Gibbs free energy of reaction can be obtained directly from the enthalpy 
of reaction: 


AGrxn(f,p) = a G rxn (T 0 ,p) _ r T Aff rxn (T» 

RT RT 0 J To RT ' 2 


(12.14) 


In order to utilize this equation, we need to know how the enthalpy of reaction varies with temperature. 


One approximation is to assume that the enthalpy of reaction is approximately independent of tem- 
perature. We can now perform the above integration with respect to T' to obtain: 

AG rxn (T, p) A G rxn (T 0 ,p) AiT rxn (T 0 ,p) T 0 \ ( /10 1CN 

RT RT RT 0 ( 1_ tJ + "' (12 ' 15) 

where AiT rxn (Tb,p) is the enthalpy of reaction at a temperature Tq. Note that this relation relies 
on the enthalpy of reaction being approximately constant in the range of temperatures from Tq to 
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T. This condition, however, can be relaxed by incorporating knowledge of the heat capacities of the 
components involved in the reaction. 


With this expression for AG rxn? the corresponding expression for the temperature dependence of the 
equilibrium constant is: 


K eq (T) 


I<eq(To) exp 


AiJ rX n(T 0 ,p) 

RTo 



(12.16) 


For an exothermic reaction, where AH Ixn (To, p) < 0, an increase in temperature will lead to a 
decrease in the equilibrium constant, thereby shifting the reaction to the reactants side. For an en- 
dothermic reaction, where AH Txn (To,p) > 0, an increase in temperature will lead to an increase in 
the equilibrium constant, thereby shifting the reaction to the products side. 
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12.6 Liquid phase reactions 


For reactions in the liquid phase, the calculations proceed in a manner similar to those for the gas 
phase. Often for liquids, it is more convenient to express the chemical potential in terms of activity 
coefficients 


Ta = V° a ( T ,P) + RT In x a -f, a 


(12.17) 


where /r° (T, p) is the chemical potential of pure species a in the liquid state at temperature T and 
pressure p. If we substitute this expression for the chemical potential into Eq. (12.7), we find 

u a p a = o 

a 

^2 Va + RT ln x ocloc) = 0 
a 

RT V a ln Xa'Ja = - yy V a P° a 

a a 

ln]T( = --Ey v a p° a 


T( x a7a) i/ “ = exp 


AGrxn^p) 

RT 


where AG rxn = u a^a(T,p) is the Gibbs free energy of reaction. 
Defining the equilibrium constant K eq as ln K eq = — AG rxn / (RT), we find 

IF^r* = K eq (T, P ) 


(12.18) 


(12.19) 


This equation is also known as the law of mass action. 
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